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ABSTRACT 

We use a simple optical/infrared (IR) photometric selection of high-redshift QSOs that identifies 
a Lyman Break in the optical photometry and requires a red IR color to distinguish QSOs from 
common interlopers. The search yields 100 z ~ 3 (J7-dropout) QSO candidates with 19 < < 22 
over 11.7 deg^ in the ELAIS-Nl (ENl) and ELAIS-N2 (EN2) fields of the Spitzer Wide-area Infrared 
Extragalactic (SWIRE) Legacy Survey. The z 3 selection is reliable, with spectroscopic follow-up of 
10 candidates confirming they are all QSOs at 2.83 < z < 3.44. We find that our z ^ A (g'-dropout) 
sample suffers from both unreliability and incompleteness but present 7 previously unidentified QSOs 
at 3.50 < z < 3.89. Detailed simulations show our z ~ 3 completeness to be ^ 80 — 90% from 
3.0 < z < 3.5, significantly better than the ~ 30 — 80% completeness of the SDSS at these redshifts. 
The resulting luminosity function extends two magnitudes fainter than SDSS and has a faint end 
slope of /3 = —1.42 ± 0.15, consistent with values measured at lower redshift. Therefore, we see no 
evidence for evolution of the faint end slope of the QSO luminosity function. Including the SDSS 
QSO sample, we have now directly measured the space density of QSOs responsible for ~ 70% of the 
QSO UV luminosity density at z ~ 3. We derive a maximum rate of HI photoionization from QSOs 
at z ~ 3.2, r — 4.8 x 10~^^ s~^, about half of the total rate inferred through studies of the Lya forest. 
Therefore, star-forming galaxies and QSOs must contribute comparably to the photoionization of HI 
in the intergalactic medium at z ~ 3. 

Subject headings: quasars — general, quasars — luminosity function: intergalactic medium 



1. INTRODUCTION 

The QSO luminosity function (QLF) is an observable 
constraint on models of galaxy formation and the corre- 
sponding growth of super-mass i ve black holes (SMBHs) 
/eg. Small & Blandford) [l99l iKauffmann fc Haehnell] 
I2000t iHaiman fc Menoul |2000D . These models are use- 
ful in interpreting observed phenomena, such as the re- 
lation bet ween a galaxy's black hole mass and bulge lu- 
minosity f Kormendv fc Richstonc 1995: Magorrian et al.l 
[1998.) ■ as well as inferring specifics of QSO activity, such 
as initial black hole mass functions, QSO light curves, 
and accretion rates (see eg. lHopkins et al.ll2006f ). 

In addition to galaxy formation models, the QLF can 
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be used to derive the QSOs' contribution to HI and 
Hell reionization. QSOs are respo nsible for H ell reion- 
ization at z 3 I Jakobsen et all IT9 94: Rci mers et al.l 
Il997t iHogan et all Il997t ISokasian et al.. ,2002 ) and are 
presumed to have a neglible contributio n to the HI reion- 
ization at z ~ 6 (jMadau et al.n i999V However, both 
of these claims require assumptions about the faint end 
slope of the QLF, as this has not been measured accu- 
rately at z > 2. 

The first QSO luminosity functions demonstrated 
a rapid incre ase i n sp ac e densities towa, r d hig her 
redshift (Schi^ [l968l : iSchmidt fc Greei] [l983l) . 
Deeper surveys, which primarily identified QSOs 
by their "UV-excess" found that the faint end 
of the QLF was shallower than the bright end 
(Bovlc et al.' 1988; Hcislcr & Ostriker|ll988|; iKoo fc Kro^ 
J^8.8; Hartwick fc S chade 1990). Recent large surveys, 
most notably the 2dF Quasar Reds hift Survey (2QZ, 
iBovle et al.H2000l : ICroom et alll2004D have found thou- 
sands of z < 2.5 QSOs. With these large samples, the 
QSOs have been placed into smaller bins in both lumi- 
nosity and redshift, accurately constraining the shape 
and evolution of the QLF. The data are typical ly fit to 
a broken power-law (|Bovle et 
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with the break at L* and a bright end slope, a, steeper 
than the faint end slope, /3. The QLF evolution with red- 
shift is consistent with Pure Luminosity Evolution (PLE) 

QSOs were most easily identified as point sources with U ~ B 
colors bluer than most stars. 
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be parametrized by a shift in the luminosity of the 
break, L*{z), without any change in its shape. Re- 
cently, the 2dF-SDS S LRG and QSO Survey (2SLAQ, 
IRichards et"alll2005[ ) has extended the "UV-excess" QSO 
search one magnitude fainter to more accurately measure 
the faint end slope of the QLF out to z ~ 2.1. These 
deeper data, when comb ined with the brig ht end slope 
from the 2QZ and 6QZ (ICroom et gd] 120041 ) . fit a faint- 
end slope, (3 — —1.45, and bright-end slope, a = —3.31, 
and demonstrates a ~40-fold increase in L* from z ~ 
to z = 2. 

At high redshift, various surveys have been con- 
ducted through gris m searches for UV emission lines 
(jSchmidt et al.lll995l ). spectroscopic follow-up of point 
sources wi th optical colo rs away from the s tella, r main 
sequence (I Warren et al.l [l994; Kcnncfick e t al.l Il995t 
iFan et ani2001l). X-ray ('Hasingcr ct al. 2005), searches 
for radio-loud (iDunlop fc Peacock 1990) or infrared lu- 
minous QSOs (iBrown et all (20061) . These surveys all 
show a precipitous decrease in space densities at z > 3. 
Early data from the Sloan Digital Sky Survey (SDSS, 
lYork et al.|[200G[ ). iFan et all (|2001i ) show a factor of six 
decrease in M < —25.5 QSOs from z = 3.5 to z = 6.0. 
Recent results from SDSS suggest that the bright end 
slope is no longer constant at z > 3. Rather, it is getting 
shallower toward higher redshift (Richards et al. 2006). 
Unfortunately, these high redshift surveys are shallow 
{i < 20) and can only measure the bright end of the 
QLF at z > 3. Therefore, little can be said about the 
shape of the faint end of the high-z QLF, or its inte- 
grated properties (eg. contribution to intergalactic HI, 
Hell ionizing radiation or black hole growth) without a 
c ensus of fain ter QSOs at z > 3. 

iHtuit et aO |2004) (hereafter, H04) searched for faint 
AGN at z ^ 3 with deep Keck spectroscopy over 0.43 
deg^ and found 11 QSOs. Though limited by small num- 
bers, the fitted faint-end slope, /3 = —1.24 ± 0.07, is 
substantially shallower than low redshift measurements. 

X-ray selected samples suggest a rather modest evo- 
lution in AGN space density at high-z (jBarger et al.l 
[2005h . or even a luminosity dependent density evolu- 
tion (LDDE) where lower luminosity AGN peak in num- 
ber density at lower redshifts than mor e luminous QSOs 
(|Ueda et al.ll2003t iHasinger et al] |2005). 

Given the importance of the QSO luminosity function 
in constraining models of galaxy and black hole forma- 
tion as well as the contribution of AGN to the ionizing 
background, these initial indications that the QLF shape 
is evolving at z > 3 warrant deeper surveys to better con- 
strain the high redshift QLF. 

The SWIRE Legacy Survey (|Lonsdale et al.ll2003f l. a 
wide-area infrared survey with deep ground-based opti- 
cal data, is optimal for searches of faint QSOs at high 
redshift as it is deep enough to detect sub-L* QSOs at 
z < 4 and covers sufhcient area to detect large numbers 
of them. In Section [3] we outline our method for creating 
a new QSO template spanning far-UV to mid-IR wave- 
lengths. In Section IH we present a simple optical/IR 
color selection for QSOs at z > 2.8 and identify areas of 
possible contaminations or incompleteness. Our selection 
results are given in Section[5l In Scction[6l the reliability 



of the sample is assessed through spectroscopic follow-up 
and analysis of the infrared colors. In Section[7]we deter- 
mine, through simulations and comparisons with known 
samples, our sample completeness as a function of red- 
shift. In Sections [8] and [9] we present our measurement 
of the QLF at z ~ 3 and compare to previous studies. In 
Section [TOl the QSO contribution to photoionization of 
HI in the IGM is computed and compared to measure- 
ments of the total photoionization rate. 

Although many early studies of high redshift QLFs 
use an Einstcin-DeSitter cosmology with Hq — 50km 
Mpc^^, throughout this paper we choose to use a 
more recent cosmology with Hq = 70 km s~^ Mpc""'^, 
fl„i = 0.3, r^A = 0.7 and correct other measurements ac- 
cordingly. All optical magnitudes are Vega magnitudes 
unless stated otherwise. 

2. OBSERVATIONS 

The SWIRE survey covers 49 deg^ over six fields at 
high galactic latitude with minimum galactic cirrus emis- 
sion. Most of this area has now been imaged in mul- 
tiple optical filters to depths r' < 24.5. Our analysis 
is conducted within the first two fields for which both 
optical and infrared catalogs were available, ELAIS-Nl 
(16hllm+55°00') and ELAIS-N2 (16h37-h41°02'). 

2.1. Spitzer Infrared Data 

SWIRE is an IR imaging survey w ith all four bands 
on the Infrared Array Camera (IRAC. lFazio et al.ll2004l ) 
and all thre e bands on t he Multiband Imaging Photome- 
ter fMIPS iRieke et al.l 12004) aboard the Spitzer Space 
Telescope (jWerner et al.ll2004f ). The IR filter character- 
istics and SWIRE depths are summarized in Table [TJ 
The ENl IRAC campaign was undertaken January 14- 
20, 2004 and MIPS January 21-28, 2004. MIPS went 
into standby mode on January 25, 2004, resulting in lost 
AORs which were reobserved July 29, 2004. The EN2 
IRAC campaign was observed July 05-06, 2004 and MIPS 
July 08-11, 2004. 

The IRAC and MIPS pho tometry were m easured with 
the SExtractor program ( Berlin fc Arnouts 1996) within 
1.9" and 15" diame ter apertures for I RAC and MIPS 
24^m, respectively (jSurace et al.l l2005l ) . Aperture cor- 
rections were derived from measurements of composite 
point spread functions from bright stars in the SWIRE 
fields. 

2.2. Optical Data 

The ENl and EN2 field s were imaged as par t of the 
Wide Field Survey fWFS. 'McM ahon et al.ll200l ". Im- 
ages were taken with the Wide Field Camera (WFC) 
on the 2.5-meter Isaac Newton Telescope (INT). Both 
ENl and EN2 have been observed with the U,g'y,i',Z 
filters over 9 deg^ each, with 600 second exposures at 
each pointing in each filter. A fraction of the fields 
~ 30% were not observed on the same night in every 
filter. The filter characteristics and depths are sum- 
marized in Table The median seeing is ~1.1" and 
never worse than 1.6". The optical coverage overlaps 
the Spitzer IR dat a by 7.45 and 4.29 deg^ in ENl and 
EN2, respectively ijSurace et al.l l2005l ). Data process- 
ing was done by the Cambridge Astronomical Survey 

http://www.ast.cam.ac.uk/ wfcsur/ 
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Unit CASU) and is out li ned in I Irwin fc Lewis! (|2001h and 
iGonzalez-Solares et afl (|2005| 1. 

Photometry was measured with the CASU software, 
requiring a source to have five continguous pixels 1.5a 
above the background pixel noise. Detection and pho- 
tometry were performed in each band separately and 
then matched between bands. Fluxes were measured 
within 2.3" (7 pixels) diameter apertures. Given the 
typical seeing of ~ 1.1", these apertures contain 80-90% 
of the total flux. Aperture corrections for each image 
were derived from bright stars within that image. To- 
tal and isophotal magnitudes were computed as well but 
only aperture photometry is used in this analysis as our 
objects are point sources by definition. Limiting mag- 
nitudes (5cr) for non-detections were computed from the 
pixel-to-pixel noise of the corresponding image. The im- 
ages were not interpolated before photometry was per- 
formed and therefore do not suffer from correlated noise 
from projection procedures. 

3. QSO TEMPLATE 

Our optical QSO selection (defined in Section |4T|) uses 
three filters which typically span the rest-frame wave- 
lengths 700 A < A < 2000 A. For the infrared selection 
(see Section |4?2)) we use the same two filters (IRACl & 
IRAC2) for all targeted QSO redshifts, so the photom- 
etry samples a large range in the rest-frame optical to 
near-IR wavelengths, 5000 A < A < 2/im. In order to 
define the expected optical/mid-IR colors of QSOs at 
2.7 < z < 5.0, we have created a QSO template which 
spans many decades in wavelength from the far-UV to 
the mid-IR. 

3.L UV Template 

Several composite QSO spectra have been created from 
large QSO surveys, and agree well with each other ex- 
cept for minor variations due to selection effects in opti- 
cal color and luminosity. IVanden Berk et all ()2001f ) cre- 
ated a composite spectrum of 2200 QSOs from the SDSS 
spanning the wavelengths between 800 A < A < 8555 A. 
In order to obtain such a broad coverage in rest-frame 
wavelength, QSOs from a wide redshift range (0.044 

< z < 4.789) were used. Since the UV spectrum is 
composed from only high-redshift QSOs, the mean con- 
tinuum shortward of Lya is artificially decreased by the 
large number of Lyman line absorbers at high redshift. 
Because the level of absorption is redshift dependent, 
we need a template which reflects the intrinsic SED, to 
which we can then apply a redshift dependent model of 
t he HI absorp t ion. 

iTelfer et al.l ()2002l ) have created a high S/N composite 
UV spectrum (300 A < A < 1275 A) using HST UV 
spectra of 77 radio quiet QSOs. Because these data were 
taken in the observed UV rather than the optical, the 
QSOs contributing to the critical wavelength range 700 

< A < 1216A are at lower redshift (< z >= 1.42) than 
the SDSS QSOs (whose spectra cover the same rest-frame 
wavelengths at z > 2.4). Therefore, there is much less 
Lyman line and continuum absorption in this composite. 
Furthermore, Telfer et al. (2002) corrected for absorbers 
with column densities Nhi < lO^^cm^ and statistically 
corrected for lower column density absorbers, resulting in 
a good template of the intrinsic UV spectrum of QSOs. 
For our UV template we have used the composite from 



ITelfer et al.l (l200l ) for 300 A < A < 1250 A and the SDSS 
composite for 2000 A < A < 8555 A. The mean of the 
two composites is used in the overlapping regions 1250 
A< A <2000 A, w here they agree well with each other 
(jTelfer et al.l[200l . 

3.2. Optical to Mid-Infrared Template 

The SDSS composite from IVanden Berk et all (|2001l ) 
covers the optical to A < 8555 A. However, at A > 
5000 A the composite is produced by low-z, low lumi- 
nosity AGN and suffers significant contamination from 
host galaxy stellar light. Indeed, a comparison with 
the median broa dband SED of luminous QSOs from 
lElvis et "all (| 19941 ) , for which stellar contribution to lumi- 
nosity should be small, shows a much redder continuum 
slope in the SDSS composite at A > 5000 A. To rehably 
select QSOs at 2.5 < z < 6.0 with our mid-IR photome- 
try, a new template is needed which covers the rest-frame 
5000 A < A < 2/im and resembles the high luminosity 
(Afi450 < —23) QSOs that we are targeting. 

We matched the SDSS photometry of all spectro- 
sco pically confirmed QSO s in the SDSS Data Release 
3 (|Schneider et al.l 120051 ) with the near-IR photome- 
try from the T wo Micron All Sky Survey (2MASS, 
ISkrutskie et al.l 2006) point-source catal og, including the 
2MA SS Deep Lockman Hole Catalog (|Beichman et al.l 
[200l . which goes ~ 1 magnitude deeper than the typical 
2MASS all-sky depths. There were 8642 2MASS objects 
within 2" of the SDSS positions. To extend the template 
to longer wavelengths, we also matched the Spitzer IR 
photometry of the 241 SDSS QSOs in the SWIRE fields 
(Lockman, ENl, & EN2). 

We chose to normalize the photometry at rest frame 
Anorrn = 2400 A for three reasons. First, the Xnorm is 
easily measured with our ground based optical/near-IR 
data for a broad range of redshifts, without using ex- 
trapolations and invoking assumptions about the SED 
spectral slope. Second, the UV flux is dominated by the 
AGN rather than the host galaxy stellar light. Third, this 
wavelength is far from any major emission lines which 
would significantly affect the photometry (eg, Mgll or 
CHI]). Only QSOs with detections in adjacent filters sur- 
rounding Xrest = 2400 A were used to minimize the inter- 
polation of the flux to the normalization wavelength. Af- 
ter an F,,(2400 A) was calculated, only those QSOs with 
-^^2400 < —24 were used to ensure the flux was dominated 
by the AGN. This also corresponds well with the mini- 
mum luminosity of QSOs in our search (M1450 = —23.5) 
and is ~1 mag brighter than the classical QSO/Seyfert 
demarcation. This cut in absolute magnitude left 3378 
QSOs (39% of SDSS/2MASS matches) with which we 
made the template. 

The template was made by averaging all of the flux 
measurements (SDSS/2MASS/SWIRE) in 1000 wave- 
length bins spaced logarithmically between 753 A and 
11 /im. The result is plotted in gre en in Figure [ 1] and 
trac es well the combined temp late of ITelfer et al.l ()2002f ) 
and IVanden Berk et al.l ()2001h in the UV. As expected, 
the relative UV to o ptical ratio of our templat e is signif- 
icantly less than the IVanden Berk et al.l ()2001h template 
because our sample is subject to less contamination from 
the host galaxies. 

Broadband photometry works well in recreating 
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Fig. 1. — The new optical/IR QSO template combined with 
ITelfer et al.l II2002I ) template (i)lack). The broadband template at 
A < 8OOOA has been replaced by the corrected Telfer / Vanden Berk 
template but is plotted in green to demonstrate the accuracy of 
the sp ectral slope. Also plotted for comparison arc thelE lyis et al.l 
1)1991 ') template (blue) and the Vandon Berk et al. (206T) SDSS 
composite spectrum (red). All templates are normalized at A = 
2900A. 

smooth, continuous SEDs, but sharp features such as 
emission hncs and continuum breaks are convolved with 
the filters through which they are observed, effectively 
broadening the features. Dcconvolution of this broad- 
ening is difficult since our photometry comes from 13 
different filters (ugrizJHK,lRACl-4,MlPS24). The only 
sharp feature at A > 5500 A is the combined emission 
lines of Ha (6563 A) & [Nil] (6548 A and 6592 A). We 
therefore use the broadband photometry to get the shape 
of the underlying continuum at A > 5100 A, subtract the 
areas affected by Ha & [NH] , and interpolate the contin- 
uum along this regi on. We th e n add the Ha-|-[NII] profile 
from the Vanden Berk et al.l (|2001h template (scaled to 
the 2400 A continuum value to preserve the relative line 
strengths). 

The resulting infrared QSO template is shown in Fig- 
ure [1] and tabulated in Table [3l An important fea- 
ture of the template is the minimum at A ~ l/im and 
the monotonic rise (in lyf^) toward longer wavelengths. 
This is produced b y hot dust at va rying distances from 
the central engine (jBarvaini i llQSl . The minimum at 
A ^ 1/xm has b een shown to be prevalent in QSO SEDs 
(jSanders et al.lil989 ) and i s attributed to s ublimation of 
dust grains at T > 1500K (|Barvaimslll987l ). 

Another important aspect of this new template 
is the increased e qui valent width (W) of Ha-|-[NII]. 
IVanden Berk et ahl (|2001h derived a combined 
W(HaH-[NH]) = 197 A from their composite spec- 
trum. Since the relative stellar contribution has been 
removed in our template, we now have W(Ha-|-[NII]) 
^ 340 A, more than 70% higher. This proves to be 
important at z ~ 4 when Ha redshifts into the mid-IR 
and significantly affects the IRAC colors. 

Finally, we point out that o ur broadband tem plate 
is significantly redder than the I Elvis et "ahl () 19941 ) tem- 
plate, with the ratio of far-UV to near-IR nearly a fac- 
tor of two larger in the latter. This is not surprising 
since their sample was selected in the soft (0.3-2.0 keV) 
and "ultrasoft" (0.1-0.3 keV) X-ray bands. Also, their 
sample is composed of more luminous QSOs, which are 
known to have bluer UV -optical colors than fainter QSOs 
(jRichards et al.ll2006a|) . 



4. HIGH REDSHIFT QSO SELECTION 

In this section we define z > 3 QSO selection crite- 
ria requiring only three optical bands and the two most 
sensitive imaging bands on the Spitzer Space Telescope 
(IRACl & IRAC2) and assess its eflacacy. The method 
consists of an optical color selection to identify a Lyman 
Break in the rest-frame UV, thereby isolating QSOs to a 
narrow redshift range (ie. z ^ 3 for U-band dropouts). 
In addition, we also require a red Spitzer IR color to 
eliminate typical contaminants in Lyman-Break Galaxy 
surveys (stars and low-z galaxies). 

4.1. Optical Selection 

The space density of neutral Hydr ogen (HI) ab- 
sorbers increases rap idly with redshift tee chtoldl 119941 : 
IWevmann et al.lll998l ). There are hundreds of absorbers 
with Nhi > lO^^cm"^ along any hne-of-sight (LOS) 
to galaxies with z > 2, resulting in Lyman line ab- 
sorption, or "blanketing" , of the source's continuum at 
Xrest < 1216 A. In addition, Lyman continuum absorp- 
tion by the absorbers significantly decreases the observed 
flux at Ar-est < 912 A. At z > 2.5, these features red- 
shift into the bluest filter accessible to ground-based tele- 
scopes (U band) and can be identified by a flux decre- 
ment at wavelengths shortward of an otherwise blue con- 
tinuum. This method, known as the Lyman Break tech- 
nique, has been used ext ensively to search for both QSOs 
and galaxies at z > 3 (llrwin et all Il991l : ISteidel et al.l 



fT996t iMadau et al.lfT99^. 

We use the lMadaul (|1995f ) prescription for the HI opac- 
ity evolution to determine average QSO colors as a func- 
tion of redshift. Figure [2] shows the U — g' , g' — r' color 
track of our our QSO template as it is redshifted and the 
HI absorption is applied. At z = 2.3 the Lyman Break 
redshifts into the U filter, resulting in redder U—g' colors, 
thereby causing the QSO to move away from the locus 
of blue, low-redshift QSOs. At z > 2.8 the QSO moves 
away from the optical color space of main sequence stars 
and into color space unoccupied by typical stars or galax- 
ies. We can therefore search for QSOs in this color space 
while minimizing contamination. 

Our color selection for both the z ~ 3 and z 4 sam- 
ple, defined in Tableland shown in Figures [2] and [3] were 
defined to select QSOs with <, 2a deviation in spectral 
slope (see Section [7.3.1|) . As we expect to eliminate stars 
and low-z galaxies from our sample with an additional 
IR selection, we are able to use a color cut which is much 
broader (and closer to the stellar locus) than optical-only 
surveys such as SDSS. 

In addition to these color criteria, a candidate QSO 
is required to be unresolved in images taken through 
the two redder filters (eg. g' and r' for z ~ 3 selec- 
tion). This minimizes contamination from low-z galaxies. 
Our sample is limited to objects brighter than r' < 22, 
corresponding to a signal-to-noise ratio SNR > 20 in 
both g' and r' . Therefore, we have sufficient sensitiv- 
ity to detect significant deviations from a point source 
(eg. FWHM > 0.5"). 

4.2. Infrared Selection 

As seen in Figures [5] and [31 we expect some contami- 
nation from stars in both our U -dropout and g'-dropout 
samples. Fortunately, the mid IR SEDs of stars are very 
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Fig. 2. — U — g' , g' — r' color-color diagram showing the 2: ~ 3 
QSO selection. The solid curve is the color track of our QSO 
template with IGM absorption applied. The black filled circles 
denote the locations in redshift increments of 0.1 from 2.9 < z < 
3.4. The shaded regions denote the color space spanned by QSOs 
with ±2(T deviations in the spectral slope (see Table [7j. The lines 
are color tracks of various galaxy templates from < z < 2: Ell 
(age=2Gyr, dotted), Sa (dashed), Sc(dot-dashed), and Sd(multiple 
dot-d ashed) taken from the GRASIL library of models l |Silva et al.l 
119981) . with the boxes corresponding to their res pective colors at 
z = 0. The black asterisks are stars from the iGunn &: Strvkeil 
l(T983l ) catalog. 
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Fig. 3. — g' — r', r' — i' color-color diagram showing the 2 ~ 4 
QSO selection. The symbols and lines are the same as in Figure [2] 

blue in color as they lie on the Rayleigh- Jeans side of the 
blackbody spectrum, whereas QSO SEDs are red, rising 
(in vfu^ towards longer wavelengths as seen in Figure [TJ 
Many groups have p roven Spitzer's ab i lity to select AGN 
with IRAQ colors (Lacy ct al. '2004 iStern et all 120051 : 
iHatziminaoglou e~al., ,2005. ). In Figure [H we plot the 
IRAC colors of all point sources with r' < 22. There are 
two clear loci of points. The objects with blue IR colors 
(lower left) are stars an d those wi t h red IR c olors (upper 
right) are AGN. Both L acv et al.l (|2004[ ) and lStern et all 
(|200a) use all four IRAC bands to select AGN, but Figure 
111 shows that the IRAC1-IRAC2 is a robust discrimina- 
tor of AGN from stars when a point-source criterion is 
also used. Also, we seek an IR selection using only IRAC 
channels 1 & 2 since they are ~ 2 magnitudes more sen- 
sitive than channels 3 & 4 given the same exposure time 
(see Table [T]). Therefore, we apply the IR color cut of 

[3.6] - [4.5] > -0.15 {AB) (2) 

in addition to our optical selection to remove stars 
from our sample. Here [3.6] & [4.5] are AB magni- 
tudes at 3.6 & 4.5 ^m respectively, corresponding to 
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Fig. 4. — IRAC colors of poi nt sources in ELA IS-Nl with r' < 
22. The dash-dotted line is the IStern et all (|2005| ) AGN selection 
and the dashed line is our selection [3.6] — [4.5] > —0.15 (AB) for 
demarcation of AGN from stars. 
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Fig. 5. — [3.6] --[4.5] color vs. redshift. The line types are the 
same as in Figure [2] Filled circles are SDSS QSOs in the SWIRE 
fields. The lighter, long dashed line is our IRAC color selection, 
[3.6] - [4.5] > -0.15(AB). 

U{A.5nIn)/f^{3.6^im) > 0.87. 

Finally, we search only sources with f^{A.5fim) > 
10/iJy (~ 7a) so that Poisson errors in the flux mea- 
surement are minimized, thereby decreasing the risk of 
contamination from sources with large errors in IRAC 
colors. We will see in Section FOI that this matches well 
our r' < 22 optical cut. 

In addition to stellar contamination, some low redshift 
galaxies may also meet our optical color criteria. This 
is due to the Balmer- or 4000A-Break being mistaken 
for a [/- or g'-dropout at z < 0.1 and 0.1 < z < 0.5, 
respectively. Our point-source criterion will remove the 
lowest redshift galaxies as they will have large angular 
diameters, but our optical images may not resolve the 
most compact galaxies at 0.3 < z < 0.5. Therefore, we 
expect to see some contamination from galaxies within 
this redshift range in the g'-dropout sample. In Figure 
m the [3.6] — [4.5] color of our QSO and galaxy templates 
are plotted versus redshift, as well as the IR colors of 
SDSS QSOs in the SWIRE fields. Unfortunately, many 
galaxies which may contaminate our g'-dropout optical 
selection (ie. 0.3 < z < 0.5) are expected to meet our 
IR color selection. This is due primarily to the presence 
of strong 3.3 lim PAH em ission in star-forming galaxies 
(jimanishi fc Dudley! l2000f ). Therefore, we expect some 
contamination from z ~ 0.4 galaxies within our z 4 
QSO sample. 

In Figure m we can see that both our QSO template, 
and all but three of the SDSS QSOs remain above our 
IRAC color cut for < z < 3.9, including the 18 of 
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19 SDSS QSOs in the targeted redshift range for U- 
dropouts. However, the QSO template dips close to our 
IR color cut between 3.8 < z < 4.8, as do four of the five 
SDSS QSOs in this redshift range. This is partly due to 
the change to bluer slope at A < I.S/ito, but is caused 
primarily by Ha redshifting into IRACl. At these red- 
shifts the Ha equivalent width is W(Ha)^ .15/im and 
the IRACl filter width is 0.7/xm. Therefore, the Ha flux 
causes the [3.6] — [4.5] color to decrease by 0.2 mags. As 
a result, the (^'-dropout selection may suffer from signif- 
cant incompleteness between 3.9 < z < 4.5 (in addition 
to the mireliability discussed above). 

4.3. Selection Summary 

• Point source brighter than r' < 22 {i' < 21.5 for 
z ~ 4 sample) to ensure proper morphological char- 
acterization and bright enough to determine a large 
U - g' {g' - r') limit. 

• IRAC2 > 10/iJy to ensure high SNR needed for 
accurate mid-IR colors. 

• Optical colors which identify a Lyman Break in the 
continuum (see Table H]). 

• Red mid-IR color ([3.6] — [4.5] > —0.15) to remove 
interlopers. 

5. RESULTS 

We performed our search in the ENl and EN2 SWIRE 
fields, covering 11.7 deg^ with both optical and IR cover- 
age. We found 100 z ~ 3 and 26 z ~ 4 QSO candidates 
which meet both our optical and IR criteria. The z ~ 3 
candidates in ENl and EN2 and their optical/infrared 
photometric data are listed in Table [5] The optical col- 
ors of the z '--^ 3 sample are plotted in Figure [6] along 
with the colors of all other point sources with r' < 22. 
We will show that the z ^ 4 sample suffers from signif- 
cant contamination and therefore do not list them. In 
Section [6l we assess the reliability of our sample through 
spectroscopic follow-up. In Section [71 we assess the com- 
pleteness through Monte-Carlo simulations and derive ef- 
fective volumes for use in computing the luminosity func- 
tion. 

6. RELIABILITY 
6.0.1. Spectroscopic Follow-Up 

Optical spectra were obtained for 10 z ~ 3 and 10 z ~ 4 
QSO candidates. Thirteen spectra for faint candidates 
(6 [/-dropouts and 7 g'-dropouts) were obtaine d with the 
Low- Resolution Imaging Spectrometer (LRIS lOke et aP 
Il995f l on the Keck I Telescope during the nights of 03-04 
March, 2005. The [/-dropout sample, with the expected 
redshift of z ^ 3, have the most prominent emission lines 
at A < 75 00 A. Therefore, on ly the blue channel (see Ap- 
pendix in lSteidel et al.ll200^ was used for these sources 
with a 1.5" wide longslit and a 300 1/mm grism blazed 
at 5000 A, giving a resolution of 1.43 A pixel ^ over 
3300 < A < 7650 A. For the g'-dropout sample, in order 
to detect the CIV (1549 A) fine, we used a dichroic at 
6800 A and used both the blue and red LRIS channel. We 
used the same grism on the blue side and the 400 1/mm 
grating blazed at 8500 A on the red side giving a blue 
side resolution of 1.43 A pixer^ from 3300 < A < 6800 





J 


/ 



























-0.5 0.0 0.5 1.0 1.5 2.0 2.5 
g'-r' (Vega) 



Fig. 6. — The U — g', g' —r' optical colors of objects with r' < 22 
and classified as point sources. The solid curve is the color track of 
our QSO template with IGM absorption applied, and the square 
denotes the z = point. The green contours denote the density 
(in color space) of objects categorized as stars by their blue IRAC 
colors ([3.6] — [4.5] < —0.15). Red crosses are point sources with red 
IRAC colors ([3.6] — [4.5] > —0.15). The blue arrows and circles 
also have red IRAC colors, but match the optical color criteria 
for 2 ~ 3 QSOs. The arrows denote upper limits where there 
is no detection in U. Spectroscopically confirmed candidates are 
circled. The photometry of two spectroscopically confirmed QSOs 
were revised and are therefore slightly out of our optical selection 
window but are still displayed here. They are both QSOs at z ~ 3 
and their photometry is included at the end of Table |5] 

A and 1.86 A pixel"! from 7000 < A < 8500 A. Total ex- 
posure times ranged from 5 to 15 minutes. Both nights 
were photometric with seeing of ~ 1.2". The spectra 
were flux calibrated with observations of the standard 
star G1919B2B from Massey et al. (1988). 

Seven additional spectra (4 z ~ 3 and 3 z 4 can- 
didates) were obtained for brighter candidates with the 
COSMIC spectrograph on the 5-meter Hale Telescope 
at Palomar Observatory on the nights of 11-14 March, 
2005. A 300 1/mm grism blazed at 5500 A was used 
with a 1.5" wide longslit, giving a dispersion of 3.1 A 
pixel"! and wavelength coverage of 3800 < A < 8500 A. 
The nights were photometric with poor seeing (^ 2 — 4") 
so long exposure times of 10-60 minutes were required. 
The spectra were flux calibrat ed with observations of the 
standard star G191B2B from lMassev et all ()1988f ). 

All ten [/-dropout candidates are QSOs with redshifts 
between 2.83 < z < 3.44, the expected redshift range 
for our sample. The spectra are shown in Figure [7l ex- 
hibiting broad Lya and CIV (1549A) lines and the QSOs 
with spectroscopic confirmation are circled in Figure [6l 
The six candidates observed with Keck/LRIS were cho- 
sen simply in ascending order in declination in order to 
sample the QSOs randomly in color space. After these 
six were shown to be QSOs, we decided to target QSOs 
with extreme optical colors on the edges of our selection 
region as these are more likely to be interlopers or ob- 
jects with spurious photometry. In our subsequent run 
at Palomar, four of these QSOs with extreme colors were 
observed: one with blue g' — r' color, one with red U — g' 
color, and two with optical colors that lie close to stellar 
main sequence. All four are in fact QSOs at the expected 
redshift. 

As expected (see Section 14.21) , our g'-dropout sample 
suffers from significant contamination. Seven of ten can- 
didates are QSOs with 3.48 < z < 3.88 and arc displayed 
in Figure [5] and listed in Table [HI The three other spec- 
tra are plotted in Figure [9] Two of the contaminants 



z=3 QSO Luminosity Function 



7 




3500 4500 5500 6500 7500 

Wovelength (A) 

Fig. 7.— Keck/LRIS (top 6) and Palomar/COSMIC (bottom 4) 
spectra of z ~ 3 QSO candidates. The ID numbers from Table \S\ 
and redshifts are given and the hja and CIV lines are labeled. 
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Fig. 8.— Keck/LRIS (top 4) and Palomar/COSMIC (bottom 2) 
spectra of z ~ 4 QSO candidates. The ID numbers from Table |6] 
and redshifts are given and the Lya line is labeled. 

are galaxies at low redshift {z = 0.354,0.390) and exhibit 
strong breaks in the continuum at A 5500 A. At these 
redshifts, the interlopers' 4000 A break falls between the 
g' and r' filters. The third contaminant has a high S/N 
spectrum, but did not show any emission features, nor a 
break in the continuum as expected, so we have not iden- 
tified a redshift for this object. Therefore, the reliability 
of our z ~ 3 and z ~ 4 QSO samples are 100% (> 69% 
la) and 70l26'^: respectively. 

The spectroscopic follow-up of the g'-dropouts shows 
the expected incompleteness at the high redshift end of 
the targetted range (3.5 < z < 4.5). All 7 confirmed 
QSOs have z < 3.9. Beyond this redshift. Ha redshifts 
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Wovelength (A) 

Fig. 9.— Keck/LRIS (top 2) and Palomar/COSMIC (bottom) 
spectra of interlopers in the g'-dropout sample. The official SWIRE 
names and redshifts (if known) are labeled. 

into IRACl and may cause the [3.6] — [4.5] color to be 
bluer than our color criterion. 

6.0.2. Infrared Reliability 

In Figure[6l there are a few tens of point sources (out of 
more than 11,000) with red IRAC colors that lie within 
the stellar locus where we would not expect QSOs to lie 
(upper-right), suggesting possible contamination where 
the stellar locus crosses the color space of our QSO se- 
lection. These may be red galaxies at moderate redshift, 
highly reddened QSOs, or stars whose IRAC colors are 
just redder than our color cut. Figure [TOl shows that in- 
deed there is not a distinct bimodality of IRAC colors for 
point sources. It's clear there are two peaks in the dis- 
tribution, but there are hundreds of sources in the valley 
between the two peaks. We have also plotted the IRAC 
color distribution of the point sources which also meet 
our optical color criteria. In this histogram, nearly all of 
the sources in the valley between the two peaks have been 
removed, revealing a clear bimodality. There are two rea- 
sons for the this. Firstly, galaxies can have either red or 
blue IRAC colors depending on redshift, but should not 
he in the optical color space of z 3 QSOs. Secondly, 
though a few stars may have IRAC colors redder than 
our cut, the stars must be either very bright or very red 
to be detected in IRAC channel 2 at all. Stars that lie in 
the z ~ 3 optical color space have r' — [4.5] = 1.6 mags 
(AB). Since the IRAC channel 2 flux limit corresponds 
to 21.4 in AB magnitudes, this means we shouldn't de- 
tect the stars in channel 2 unless they are brighter than 
r' < 19.7 (Vega). Indeed, in Figure im there are two 
sources at the bright end of our sample that have blue 
optical-IR colors and may be stars. 

As mentioned in Section 14.21 IRAC color selection of 
AGN is robust when all four bands are used because the 
SEDs of most contaminants are not flat (or rising toward 
longer wavelengths) over such a broad wavelength range. 
Therefore, we can use all four IRAC fluxes, when avail- 
able, to assess the reliability of our single color AGN 
selection. In our sample, 58% (58 of 100) of our ob- 
jects are detected in all four IRAC channels and their 
IRAC colors are plotted in Figure [T^l Essentially all 
of these objects have IR colors in the expected locus of 
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Fig. 10.— The IRAC ([3.6] - [4.5]) color distribution of all point 
sources with r' < 22 in our field (dashed line). The solid line 
show the color distribution of point sources which also meet our 
optical color criteria. The shaded histogram is the distribution of 
spectroscopically confirmed QSOs. Our QSO sample is the portion 
of the solid histogram to the right of the vertical dotted line. 
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Fig. 11. — Optical-IR color distribution of the z ~ 3 QSO sample. 
The X-axis is in Vega magnitudes since our magnitude bins are 
defined in Vega magnitudes and the y-axis is in AB magnitudes as 
it is simpler to interpret. The dashed line in the top plot denotes 
the /4.5^m = lO^Jy limit of our search. The dashed line in the 
bottom plot denotes the /24^m = 250/im completeness limit with 
arrows giving upper limits from non-detections. The four QSOs 
with r' — [4.5] < 0.0 are possible interlopers. 

AGN. Therefore, it appears the contaminatfon rate is 
fow. However, as steUar contaminants are expected to 
have blue IRAC colors, we do not expect to detect the 
fainter interlopers in IRAC3 and IRAC4 and we cannot 
assess the nature of these objects. Therefore, we look 
at only the bright end of our sample where all sources 
are detected in all four bands. At r' < 20.5, 21 of 24 
(88%) sources are detected in all four IRAC bands, and 
all of these have IR colors within the expected AGN lo- 
cus. If we make the conservative assumption that the 
three non-detections (in IRAC3 & IRAC4) amongst the 
bright sample are not QSOs, then we obtain an upper 
limit to the contamination rate of <13%. 

In summary, the point sources which match our op- 
tical color criteria have a distinct bimodal distribution 
in [3.6] — [4.5] colors, suggesting minimal contamination 
from stars, except at the bright end (r' < 20). However, 
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Fig. 12.— IRAC colors of the z 3 QSO sample. The lines 
are as in Figure |4]. Red arrows are lower limits to the [5.8] — [8.0] 
color based on non-detections in IRAC3. Blue arrows are upper 
limits to the 5.8] — [8.0] color based on non-detections in IRAC4. 
Non-detections in both IRAC3 & IRAC4 are plotted along the left 
hand side. 

the four-band IRAC colors suggest that even the bright 
end of our sample has a contamination rate of less than 
13%. 

7. COMPLETENESS 
7.1. Infrared Completeness 



As mentioned in Section 14.21 we expect our IR color 
cut to include nearly all QSOs at z '--^ 3 (as seen with 
the SDSS QSOs). However, we must ensure that no 
QSOs are missing due to the /^(4.5/im) > lOfiJy. Fig- 
ure [TT] plots the optical-IR colors of all U -dropout can- 
didates as a function of the optical magnitudes. The lo- 
cu s of colors ag r ee wel l with the infrared selected QSOs 
of iBrown et al.l ()2006l ). Also plotted is the limit due 
to the IRAC2 flux cut. Except for the faintest of the 
six half-magnitude bins, there is is no incompleteness 
due to non-detections in IRAC2. In the faintest bin, 
21.5 < r' < 22.0, it is possible that we miss a few of 
the bluest r' — [4.5] QSOs, but we estimate this to be 
< 15% (based on the r'-[4.5] distribution in Figure fTTj) 
and make no attempt to correct for it. We've also veri- 
fied that increasing the counts by 15% in this bin has a 
neglible effect on the QLF fit. 

7.2. Morphological Completeness 

Our selection requires that the candidates be unre- 
solved in the optical data, where typical resolution is 
^ 1.1". We have assumed that the QSOs are signifi- 
cantly brighter than their host galaxies and, even if we 
could detect the host galaxy, galaxies at z > 3 shoul d 
have angular diameters less than 1" (jCiavaliscol l2002| ). 
To test this assumption, we've matched our catalogs to 
all of the SDSS QSOs (which do not have to meet a mor- 
phological criterion) within the ENl and EN2 fields. Of 
the 58 SDSS QSOs at high redshift (z > 1), 57 (98%) 
are categorized as point sources in our optical catalogs. 
Therefore, we do not expect any significant incomplete- 
ness due to our point-source criterion. 

7.3. Optical Completeness 

Our color selection should be strict enough to min- 
imize contamination, but broad enough to encompass 
the majority of the targeted QSOs to minimize com- 
pleteness corrections. As a first-order look at com- 
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pleteness, we compare with the SDSS QSO sample be- 
tween 3.1 < z < 3.2 w here SDSS is ^ 80% completete 
(|Richards et al.ll2006bD . SDSS uses four colors, to look 
at all point-sources a way from the color space of stars 
(jRichards et al ."2002i. Therefore it is possible that SDSS 
selects redder QSOs at z ~ 3, which would lie outside 
of our color selection. Converting SDSS to Vega mag- 
nitudes, and correcting for slight U/u' filter differences 
(0.06 mag correction), -97.5% (392/402) of SDSS QSOs 
at 3.1 < z < 3.2 would also be selected with our color 
selection. Therefore, < 3% of SDSS QSOs are redder, 
and selected with other colors in the SDSS filters. 

Significant dispersion in spectral features may cause 
some QSOs to lie outside of our selection criteria. The 
most important factors affecting the optical color (and 
therefore the completeness) are redshift, intervening high 
column-density HI absorbers, UV continuum slope, and 
emission line equivalent widths, and variability. Here we 
present our Monte Carlo simulations to assess the com- 
bined effects of these characteristics on our completeness 
as a function of magnitude and redshift. We then use 
this completeness to derive the effective volume of our 
sample as a function of apparent magnitude. 

7.3.1. Model of QSO Optical Color Distribution 

Both the UV spectral slope of the continuum a^, where 
fi, oc , and the emission line equivalent widths, W, can 
vary about the mean and cause significant dispersion in 
colors at any given redshift. The mean of the spectral 
slope and the equivalent widths are defined by our tem- 
plate spectrum and are co nsist ent with values fou nd by 
IVanden Berk et all (|2001h and iHunt etlHI (|2004D . The 
distribution of these attributes are assumed to be gaus- 
sian with standard deviations taken from Francis (1996) , 
and are listed in Table [T] 

Along the line-of-sight to any given QSO at high red- 
shift there are hundreds of intervening neutral hydrogen 
clouds (or filaments) absorbing the UV light through 
Lyman line and Lyman continuum absorption. Much 
of this absorption is caused by the rare, high column- 
density absorbers known as Lyman Limit Systems (LLSs, 
Nhi > 1.6 X lO^^cm^) and Damped Lyman-a Systems 
(DLAs; Nhi > 1 x 10^° cm^). Therefore, there is a large 
dispersion in line-of-sight HI opacity at any given red- 
shift due to the small numb e rs of these high column den- 
sity clouds. iBershadv et al.l ()1999i) show that simple ana- 
lytic expressions for the mean and scatter in line-of-sight 
opacity are insufficient and, therefore, Monte-Carlo sim- 
ulations are required to correctly represent the stochastic 
distributions of the absorbers. 

The number density distribution of HI absorbers is typ- 
ically given as a fit to a power-law with redshift 

N{z) = Noil + zy (3) 
and the column density distribution is given as 



Since absorbers of d ifferent column densities are known 
to evolve differently ('Ki m et al.ll9 97^. we have split them 
into fo ur groups with d i fferin g evo lution, similar to the 
work of lBershadv et all (|1999[ ) and lFan et all (|2001h but 
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Fig. 13. — A subsection of the simulated IGM transmission 
(e~^) curve for one simulated line-of sight. All aborption features 
are fully resolved in the simulation. 

with updated values for DLA number densities. The pa- 
rameters used are summarized in Table [51 

For each line of sight, the number of line-of-sight ab- 
sorbers is selected from a Poisson distribution with the 
expectation value set to < TV > of the population. 

rZQso 

<N >^ N{z)dz (5) 

Jo 

Their redshifts and column densities are chosen randomly 
from the distributions in Eqns. [3] and 31 respectively. Fi- 
nally, we compute Voigt profiles with natural broadening 
and doppler widths, &, from Table [51 for the first ten Ly- 
man lines for each absorber. We choose not to model the 
distribution of Doppler widths (or its evolution) as this is 
expected to be a small, second order effect. Continuum 
absorption is applied with a scattering cross-section of 
(To = 6.3 X 10^^^ cm^ at the Lyman Limit and decreas- 
ing as v"^ . An example of the Lya Forest transmission 
for one line of site at z = 3 is plotted in Figure [131 

For 300 lines-of-sight in each redshift bin of Az = 0.1, 
we compute QSO spectra with spectral slopes and emis- 
sion line equivalent widths culled from the distributions 
defined in Table [51 and apply the simulated IGM absorp- 
tion for that line of sight. 

QSOs of these luminosities and redshifts exhibi t vari- 
abilit y of order Am ~ 0.15 mags (Vandcn Berk et al.l 
I2OOI . 29% of our fields were observed at separate times 
(over two years) in the U and g' filters, and 29% of our 
fields were observed at separate times in the g' and r' 
filters. We therefore add a variable offset selected from a 
gaussian distribution with Am — 0.15 to the computed 
U — g' and g' — r' optical colors of 30% of the simulated 
lines of sight. The net effect of this variability is a slight 
~ 5% decrease in completeness, and a slightly broader 
redshift range. 

After adding in variability and photometric errors, we 
compute the optical colors for the 300 QSOs in each red- 
shift bin. This allows us to determine the selection com- 
pleteness based on color selection alone, but the imaging 
depth must still be taken into consideration. 

The optical magnitude limit of our QSO search, r' — 
22.0, is set by our infrared depths and corresponds to 
~ 20(7 detections in both g' and r' . Therefore, we do not 
expect any incompleteness due to non-detections in g' or 
r' . However, when looking for significant "[/-dropouts", 
it's important to have deeper imaging in U relative to the 
bands at longer wavelengths. Otherwise, non-detections 
will result in upper limits in flux which cannot distinguish 
between a true "dropout" or a red source which is just 
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Fig. 14. — The U-band depth map for ENl. The grey-scale 
represents U magnitude Hmits in intervals of 0.1 mags from 22.5 
to 24.0 (Vega) with most of the area between 23.0-23.5 

below the detection limit. The WFS was not designed 
for finding [/-dropouts and the U band is generally less 
sensitive due to CCD quantum efficiency (QE) and poor 
throughput in the telescope optics. Therefore, the U- 
band images in the SWIRE fields are less sensitive (to all 
classes of object) than the g' and r' images. Furthermore, 
this sensitivity varies by ~ 0.5 magnitudes from pointing 
to pointing due to changing observing conditions (seeing, 
airmass, lunar phase). Therefore, careful measurements 
of the selection completeness as a function of pointing 
must be assessed, along with dependencies on redshift 
and magnitude. 

Both ENl and EN2 are comprised of 54 22.8' x 11.4' 
optical pointings. For each CCD at each pointing, a 5a 
limiting magnitude is determined within our 2.3" diame- 
ter aperture and is used as an upper limit when no object 
is detected. 

Due to overlaps in pointings and varying observing con- 
ditions, it is nontrivial to compute the limiting magni- 
tude for a given position on the sky. For both fields 
we make "depth maps" by making mosaic images of the 
entire field with pixel values set to the ba limiting mag- 
nitude of the deepest pointing which covers that pixel. 
In this way, we can quickly compute the area over which 
we can detect an object of a given magnitude Utest by 
summing the pixels with value greater than Utest- These 
depth maps are trimmed to the area which also has 
Spitzer data. The U-band depth map for ENl is plotted 
in Figure I14| showing the non-uniformity and complex 
structure of the coverage. 

As shown in Figure [H the bluer of the two optical 
colors, {U — g'), becomes redder as the QSO goes to 
higher redshift, making non-detections in U increasingly 
likely at higher redshift. For example, a typical r' = 
21 QSO at = 2.9 will have g' ~ 21.5 and U - 21.8, 
bright enough to be detected in all bands by our survey. 
However, a QSO with the same r' magnitude aX z = 3.4 
will have g' ^ 22 and U ~ 24.5, too faint to be detected in 
our U images. Our upper limit in this case is insufficient 
in distinguishing a high-redshift QSO candidate from a 
low-redshift object with bluer U — g'. 



Fig. 15. — The 2 ~ 3 QSO completeness C{r',z) contours as a 
function of apparent r' magnitude and redshift. The contours are 
spaced at AC = 0.1 intervals. 

Given our simulated spectra and our depth maps, we 
determine an effective completeness in the following man- 
ner. For each redshift bin, we compute U — g' and g' — r' 
colors for the 100 simulated QSOs along 100 different 
lines of sight. Then we compute the percentage of these 
QSOs which would be selected by our color criteria if our 
imaging was sufficiently deep. This gives us a measure 
of our completeness based strictly on our color criteria 
alone, C color- The completeness cutoff at the high red- 
shift end (z ~ 3.5) is nearly a step-function since the 
QSO color track moves perpendicular to our color cuts 
and most of the dispersion in color is parallel to the color 
cut. At the low redshift end (z ^ 2.9), the incomplete- 
ness is predominantly due to line-of-sight variations in 
IGM. 

Of the QSOs which meet the color criteria, the fluxes 
are scaled to give the desired r' magnitude in intervals 
of Ar' = 0.25 mags and the g' — r' color is used to de- 
termine the [/-band depth required to either detect this 
QSO or derive a lower limit to the magnitude which is 
high enough to put it in the color-color selection window. 
The percentage of QSOs with colors in our selection win- 
dow which would also be selected given the U -band depth 
at that pixel value is denoted as C(r', z), plotted in Fig- 
ure [151 The effective volume of the survey can then be 
calculated as 

dV 

Veffir')=dn C{r',z)—dz (6) 

J z=0 

where Vt is the solid angle of the survey and dV/ dz is the 
differential comoving volume. The effective volumes and 
average redshifts for each half-magnitude bin are given 
in Table [H To give an idea of the scales of the incom- 
pleteness corrections, the effective volume in our faintest 
magnitude bin Ve//(r' = 21.75) is 74% of the effective 
volume in our brightest bin Veff{r' = 19.25), requiring a 
relatively small correction of 35% to the number counts. 

8. THE Z ~ 3 QSO LUMINOSITY FUNCTION 

Given the effective volume, Veff{r'), the comoving 
space density is then 

. ,7) 

Veff{r) WMn 

where N{r') is the number of QSOs in the r' bin and 
Wbin is the width of the bin in magnitudes. We then 
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convert r' to the absolute AB magnitude at 1450A since 
the r' filter covers 1450 A over our entire redshift range 
and this value is generally used for high redshift QSO 
studies. 

Mi45o = r'+r'AB{yega)-DM{z = 3.2)+2.5/o.g(l+3.2)+is:(A 

(8) 

where r' is the Vega magnitude listed in Table O 
r'j^g{Vega) = 0.15 is the Vega to AB conversion, DM = 
47.19 is the distance modulus at z=3.2, and K{Az) is 
the K-correction resulting from shifts in redshift around 
z = 3.2. As we don't have exact redshifts for most of our 
QSOs, we set K{Az) = 0, but note that this varies by 
±0.1 mags from 2.9 < z < 3.5. 

The resulting QSO luminosity function at z ~ 3.2 
is plotted in Figure 1161 along with previous surveys of 
QSOs at these redshifts. Our measurements, while using 
smaller magnitude bins {wbm — 0.5 mags) than previous 
surveys of faint high-z QSOs, have significantly reduced 
error bars. The space densities at the bright end of our 
survey match well with the faintest bins from the SDSS 
(jRichards et al.l[2006bf l and show a clear transition to a 
shallower slope at the f aint end. We use our data in com- 
bination with the .Richards et al.l (|2006b[ ) SDSS results at 
z ~ 3.25 because it is the largest sample available, and its 
completeness corrections have been carefully determined. 

We do a least-squares fit to the standard double power- 
law given in Eqn. [U converted to absolute magnitudes 



$(Ml450,z) = 



0.92 X ^Mi45o) 



(9) 

As mentioned previously, when the error bars are re- 
duced, the quasar luminosity function exhibit s curvature 
over all luminosities, not i ust at the break (jWolf et al.l 
[2003t iRichards et al.l[2005h . Therefore, since both our 
data and the SDSS data show some curvature, and the 
SDSS has much smaller error bars, a least squares fit of 
a double power-law will force the break position (M*) 
to be contained within the range of magnitudes covered 
by the SDSS. This produces a very steep faint end slope 
and is a poor fit to our faintest bins. Of the four pa- 
rameters that determine the QLF, the SDSS can only 
measure the bright-end slope with certainty at 2; ~ 3.2. 
Therefore, we also perform a fit with a fixed bright end 
slope, a = —2.85, defined by the SDSS measured rela- 
tion of a with redshift at z — 3.2. This yields a more 
realistic values of M* — —25.6, and a faint end slope, 
P = — 1.62±0.19, which is constrained by 4-5 bins fainter 
than M*. The results of the double power-law fits are 
given in Table [TOl 

8.1. Maximum Likelihood Fit 

In the above fit to the binned data, we assumed that 
all of our QSOs are at a single redshift (z = 3.2). For 
most of our objects, however, we do not have spectro- 
scopic redshifts and do not know the absolute magni- 
tudes. This can be problematic because we do not know, 
for example, if the brighter QSOs are indeed more lu- 
minous or simply at the low redshift end of our redshift 
range (the distance modulus changes by ^ 0.5 mags from 
z — 2.8 — 3.4). Furthermore, the luminosity function is 
known to change over these redshifs, which may skew the 



apparent magnitude distribution from what is expected 
at a single redshift. 

Our only known quantity is the distribution of apparent 
magnitudes of QSOs in this redshift range. Therefore, we 
have modeled the expected distribution in QSO apparent 
magnitudes to compare to the observed distribution. The 
model distributions are computed in the following way. 

1. We allow $*, M*, and faint-end slope, /3, to vary 
independent of one another. We keep the bright- 
end slope fix ed at a = —2.85 as determined by 
SDSS QSOs (IRichards et al.l[2006bh . 

2. For each set of parameters, we compute the appar- 
ent magnitude distribution in small redshift inter- 
vals. 

3. At each redshift, we apply the completeness func- 
tion as a function of apparent magnitude, C(r', z), 
computed in Section [7.3.1l 

4. Finally, we sum up the apparent magnitude dis- 
tribution in each redshift interval to determine the 
expected apparent magnitude distribution function 
over the entire redshift range for each set of QLF 
parameters. 

5. We repeat the above steps three times with differ- 
ent QLF evolution: no evolution, pure luminosity 
evolution (PLE) and pure density evolution (PDE). 

For the two evolving models, we choose the level of 

30 ) ' evolution to fit the variation in space density of 
bright QSOs seen by SDSS ('y 40% decrease from 
z - 2.8 - 3.5. IRichards et al.ll2006bh . 

In addition to the observed apparent magnitude dis- 
tri bution from ou r sampl e, we also use the QSO sample 
of IRichards et all (|2006bD from the SDSS Data Release 
3. We only use the QSOs with 2.9 < z < 3.5 and use 
the com pleteness for each QSO derived in 'Richards et al.) 
(I2006bn . 

The apparent magnitude distribution function gives 
the relative probability of finding a QSO with a given 
magnitude from the whole sample. We compute the like- 
lihood that a given parameter set adequately describes 
the data as the product of the values of the apparent 
magnitude distribution func tion for all of the Q SO r' 
magnitudes in our sample ([Marshall et al.lfl983l ). We 
then find the set of parameters which maximizes this 
likelihood. The best fit model parameters are given in 
Table [101 The maximum-likelihood fit gives a some- 
what shallower faint-end slope, P = —1.42 ± 0.15, and 
a similar location of the break, and agrees within ^ la 
with the binned QLF when we assumed all QSOs were 
at z = 3.2. The apparent magnitude distributions which 
included evolution in the QLF (both PLE and PDE) pro- 
duce best-fit parameters which were not significantly dif- 
ferent from the no evolution fits because there is little 
evolution over such a small redshift range. We choose 
to use the PLE model, with P = -1.42 ± 0.15, in the 
subsequent analysis. 

9. COMPARISON WITH PREVIOUS LUMINOSITY 
FUNCTIONS 
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Fig. 16.— The z ~ 3.2 QSO luminosity function. The 
binned data are plotted from this work (black cir c les and black 
arrows for upper limits), SDSS Richards et al.f (I2006 bj, blue 
squares), COMBO-17 (Wolf ct al. 2003, red diamonds). Hu nt et al.l 
||2p04, cyan triangles), WHO (iWarren et al.lll994 . green asterisks), 
[Fbn tanot et al. (2007, open purple circles), and Bongiorno et al. 
l2007i . orange stars). Filled symbols were used in the fitting of 
the QLF. Also plotted are the fitted QLFs using the binned data 
(dot-dashed line), the maximum likelihood fit with Pure Luminos- 
ity Evolutio n (solid line), and the QLFs of Pei (1995, dotted line), 
I Hunt et all p004 . magenta dashed line), and i Wolf et al.i 1)20031 . red 
long dashed line). 



Plotted in Figure [16] are thr ee previou sly determined 
luminosity functions at z ~ 3. iPeil (|1995f ) compiled sev- 
eral QSO samples to produce a QLF and its evolution 
between < z < 4.5. A constant QLF shape was as- 
sumed, and fitted with Pure Luminosity Evolution. Un- 
fortunately, most of the high rcdshift QSOs are very lu- 



minous, so the faint-end slope is mostly det ermi ned b' 

lam 



I 



lower redshift QSOs. As seen in FigurefTBl thell^e: 
QLF has a steeper faint end slope, (3 — —1.6 and hes 
above our QLF determination until M^^^q < —26, where 
it passes through the SDSS data points. We believe that 
this discrepancy is caus ed by the assumption of a con- 
stant QLF shape. The [Pel ()1995l ) QLF shape was de- 
termined mostly by QSOs at z < 3, and normalized to 
the bright QSOs at high-redshift. Because the bright 
end sl ope at z > 3 is getting shallower (Richards et al] 
l2006bf ). this has caused a significant over-estimate of the 
space densities of faint QSOs. In fact, all previous QLF 
estimates at high redshift that assume a shape derived 
at low redshift and normalize using bright, QSOs will 
overestimate the number of faint QSOs at z > 3. At 
z ~ 3, this results in a factor of ^ 2 overestimate, but 
will get worse at higher redshift as the bright end slope 
is measured to get even shallower at 3 < z < 5. 

The H04 QLF, though measured with a much smaller 
sample, shows a very shallow faint-end slope and a steep 
bright-end slope, resulting in a very prominent break. 
Our QLF has a significantly shallower bright-end slope 
and a steeper faint end slope. The space densities of our 
faintest QSOs are about twice that of H04 and predict 
2-3 times more QSOs between —24 < M1450 < 20. 

H04 threw out most of their AGN sample (16 of 29) be- 
cause the emission line widths were less than 2000 km/s. 
These narrow-line AGN are, on average, one magnitude 
fainter than the broad-line sample. Therefore, if these 
were included, they would have added significantly to the 
faint-end counts, resulting in a steeper faint-end slope. 
We cannot discriminate between different types of AGN 
in our sample, since we only have spectra for 10 objects 



(9 of 10 have FWHM > 2000 km/s). However, narrow- 
line AGN (ie., FWHM < 2000 km s'^) may explain at 
least part of the discrepancy in our faint-end slopes. 

Recent surveys of X-ray selected AGN have also con- 
cluded t hat the faint end s l ope gets shallower at higher 
redshift (lUeda et al.l [200l iHasineer et all 120051). The 
iHasinger et al.l (|2005f ) sample overlaps our redshift range 
but differs from ours in the same way as H04, as it is a 
soft X-ray selected typ e-I AGN sampl e . If th e difference 
between the II04 and Hasinger et al.l 1)200 5') luminosity 
functions and our QLF is attributed to increasing num- 
bers of moderately obscurred AGN at fainter re s t-fram e 
UV luminosities (as suggested by lUeda et all ([2003)), 
then we would expect to see increasing numbers of QSOs 
with redder UV spectral slopes and higher Ljji/ Lijy ra- 
tios amongst faint QSOs. Indeed, Figure [Tl] shows a 
population of r' > 21 QSOs which are redder, in both 
r' — [4.5] and r' — [24], than any at r' < 21, which may 
not be included in the II04 or IHasinger et all (j2005) LFs. 

A recent optical/X-ray se arch for faint AGN at 3.1 < 
z < 5.2 by iFontanot et al.l (j2007) finds relatively large 
space densities, requiring a steep faint-end slope, (3 — 
— 1.71 ± 0.41. Though the errors are large, this rules 
out a significantly shallower faint-end slope at z > 3 and 
agrees, within la, with our fit. In their study, only 18% 
(2/11) AGN have narrow lines. Therefore, the fraction of 
narrow-line to broad-line AGN is too small to completely 
explain the discrepancy between the shallow faint-end 
slope of H04 and the steep f a int-en d slope found in this 
study and in Fontanot et al. ("20071). 
The COMBO-17 JWoff et aU 



[200I and VVDS 



(jBongiorno et all 12007^ QLFs IBongiorno et al.l (|2007f ) 
are both consistent with our QLF, within errors. Both 
of these surveys did not require any optical color or mor- 
phological crite ria yet t hey still agree wit h our num- 
bers. Both the (jWolf et al. 2003) and (Bon giorno et all 
l2007h samples are larger than ours (192 and 130, respec- 
tively), but the QSOs are spread out over all redshifts 
(0 < z < 5) so they have 5-10 times fewer QSOs in this 
redshift range. 

10. QSO CONTRIBUTION TO HI IONIZING FLUX AT 
Z ~ 3.2 

With this new determination of the QSO Luminosity 
Function, we can simply integrate ^(L) given in Equa- 
tion [1] to determine the specific luminosity density (at 
A = 1450A) of QSOs at z - 3.2, 



e = J $(L) L dL. 



(10) 



Integrating from -30 < M1450 < -20, (43.96 < 
log(Li45o) < 47.96 ergs s~^), we derive a value for the 
specific luminosity density, 61459 = 7.3 x 10^** ergs s~^ 
Hz^^ h Mpc""^. This is comparable to the values derived 
from the iHunt et all ^OOA) QLF (7.1 x 10 ^^) an d sig- 
nificantly lower than the value from the iPeil ([1995) QLF 
(1.4 X 10^^) when correcting for different cosmologies. Al- 
though our QLF determination predi cts more integrate d 
UV flux from faint QSOs than that of'Hu nt et all (|2004D . 
our Al* is more than 1.5 magnitudes fainter. These two 
effects essentially cancel out to produce similar luminos- 
ity densities. 

We can now determine the photoionization rate 
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di^An — — an I M- 



(11) 



However, we can not assume that J{y) is the shape of the 
average QSO SED, since the higher energy photons wiU 
be rep rocessed by HI and HeH, resulting in a higher value 
for r. iHaardt fc Madaul (|1996h have modeled this repro- 
cessing in a "clumpy" IGM to determine the effect on the 
HI photoionization by QSOs. This correctly includes HI 
clouds as sources of ionizing photons, as well as sinks, and 
increases F by ~ 40% at z ~ 3. We multiply £1459 by the 
ratio of /912//1450 = 0.58 in our template, and convert 
to a proper volume emissivity a t z = 3.2 to directly com - 
pare with the eg calculated bv Haar dt fc Madaul ()1996f ) 
and how it scales with the photoionization rate, Thi and 
the ionizing intensity at the Lyman Limit J912. It should 
be stated that the scaling relations from eg to J912 and 
Thi are dependent upon the value of eg itself, as this will 
effect the ionization levels of the surrounding medium. 
However this is a secondary effect, an d our cq is within 
~ 50% of the IHaardt fc Madaul (|1996D value so we don't 
expect this to be a large effect on the scaling relations. 

We get values of Thi ~ 4.5 x 10~^^ s~^ and J912 = 
1.5 X 10~^^ ergs s ~^ cm~^ Hz~^ sr~^. As pointed out by 
iHunt et al] ()2004D , these should be taken as upper limits 
as it assumes that all ionizing photons escape from QSOs 
of all luminosities, though this is not at all clear for lower 
luminosity AGN. 

It is interesting to compare the HI photoionization rate 
from QSOs with that from star-forming galaxies (LBGs 
at these redshifts). Unfortunately, it is difficult to deter- 
mine the photoionization rate from Lyman Break galax- 
ies, since their photoionizing SEDs are difficult to di- 
rectly detect and are sensitive to parameters with large 
uncertainties: dust reddening, initial mass function, star- 
burst age, metall icity, and the esca pe fraction of ioniz- 
ing photons /p.,^(|Steidel et al.ll2001l : IShaplev et al.ll2006l : 
ISiana et al.ll2007f) . 

It is possible, however, to address this question indi- 
rectly if the total ionizing background (QSOs-f galaxies) 
is accurately determined. Several groups have made 
these measurements by measuring the mean trans- 
missi on of QSO UV flux through t h e Lyg for- 
est (McDonald fc Miralda-Escudc 2001; 'Tvtier et all 
■2004; Bohon ct al. 2005; Bohon & Hachnclt 2007; 
Becker et al.ll2007|). or by measuring the extent of the 

around high red- 
Figure [T7| shows 



proximity e ffect (ICarswe 



shift QSOs (jScott et al.l 



l_et al."l982 
2000, 2002). 



the current estimates for the total photoionization rate 
at high redshift, compared with our estimate of the 
contribution from QSOs at z ~ 3.2. The lower limit 
to F is derived by integrating the QLF between 30 < 
M1450 < 23.5, the range covered by our survey and 
SDSS. The upper limit is derived by integrating the 
QLF 3.5 magnitudes fainter and assumes a 100% es- 
cape fractio n amongst these faint QSOs as well. The 
IScott et all (12000) value has been scaled (by the value 
given in [Scott ct al. (2002 )) to our assumed cosmol- 
ogy. The ,McDonald fc Miralda-Escudel (|2001h values 
have also been scaled to the same cosmology using Eqn. 
3 in their paper. 

Although the error bars are large, all of the total 
photoionization measurements are consistent and give a 




Fig. 17. — HI photoionization rate per atom (in units of 1 X 10 ^"^ 
s~^) versus redshift. Estimates of t he total photoioniza tion rate of 



the IGM are plotted as asterisk s ITvtler et al.l 12004), diamonds 
(McDonald & Miralda-Escudc 200 1I T! and squares l lBolton et al.l 
2005.; Bolton k Ha chnclt 2 0071 ) . The shaded reg ion is the pho- 
toionization rate determined bv lScott et al.| l|2000l) ad j usted to our 
cosmology (decrea se of 31%) a s stated in S cott et al.l l|2002l) . The 
open circle is the IHunt et al.l ((2004) determination of the QSO 
contribution at 2 ~ 3. The black square encompasses the redshift 
range and plausible limits of the QSO contribution to HI pho- 
toionization from our sample. The lower and upper bounds are 
determined by integrating the QSO luminosity function to ij < 22 
(Afi450 < -23.5) and R < 25.5 (Afi450 < -20), respectively 



value of F - 1.0 X lO^^ g-i 



3. The contribu- 



tion from QSOs is less than half of this value. There- 
fore, it is likely that star-forming galaxies' contribution 
to the HI photoionization rate is comparable to that of 
QS Os at z ~ 3 . 2. Th is is consistent with measurements 
by iShull et al.l (I200I which examine the relative rates 
of Hell and HI photoionization within individual Lya 
absorbers at 2.3 < z < 2.9 to infer the spectral index 
of background radiation at each location. They con- 
clude that the spectral index varies greatly between ab- 
sorbers, with significant contribution from "soft" sources 
that may be starburst galaxies or dust-attenuated AGN. 

Furthermore, recent measurements of the escape frac- 
tion, fesc, of photoionizing radiation from Lyman Break 
Galaxies also suggest that star-formation may signifi- 
cantly contribute t o the ionizing background at z ^ 3. 
ISteidel et al.l (|2001h made a composite rest-frame ultra- 
violet spectrum of 29 LBGs and found that greater than 
50% of the photoionizing flux that is not absorbed by 
dust escapes into the IGM (ie. relative escape fraction 
fesc,rei > 0.5). Deeper spectra of 14 z ~ 3 LBGs give 
a smaller value of fesc,rei — 0.14, but this still gives an 
ionizing radiation field J900 ^ 2.6 x 10^^^ erg s~^ cm~^ 
Hz^^, nearly twice the value of our upper limit from 
QSOs. 



11. SUMMARY 

We present our method of finding high redshift, z > 
2.8, QSOs by identifying a Lyman Break in the optical 
photometry, and ensuring red mid-IR ([3.6] — [4.5]) colors 
indicative of QSOs. The use of only three optical filters 
allows a search over larger areas in the SWIRE fields as 
most of the area does not have coverage in four or more 
bands. The use of only IRACl and IRAC2 channels is 
emphasized as these two bands are a factor of seven times 
more sensitive than IRAC3 and IRAC4. 

Spectroscopic follow-up of 10 z ^ 3 ([/-dropout) candi- 
dates confirms that all 10 are QSOs between 2.83 < z < 
3.44. Spectroscopy of 10 z ^ 4 ((^'-dropout) candidates 
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confirmed 7 QSOs with 3.48 < z < 3.88, two galaxies 
at low redshift (z = 0.354, 0.390) and one unconfirmed 
redshift. We place reliability estimates on our z ~ 3 and 
z ~ 4 samples of 100% (> 69% la) and respec- 
tively. Since we have not spectroscopically confirmed all 
of our candidates, we only use the more reliable z ~ 3 
sample for determining a luminosity function. 

By using detailed models which include variations in 
number and column density of line-of-sight HI absorbers, 
UV spectral slope, emission line equivalent width, red- 
shift, observed magnitude, and photometric errors, we 
assess the completeness of the optical color selection. 
Completeness near the center of the redshift range of our 
[/-dropout selection is 85-90%. However, our complete- 
ness decreases significantly in our faintest magnitudes 
bins (~ 75%), due to the shallow depth of the [/-band 
imaging. 

We find 100 z - 3 QSO candidates with r' < 22 over 
11.7 deg^. Through our models of completeness ver- 
sus redshift, we derive effective volumes for each half- 
magnitude bin and compute the z ~ 3 QSO lumi- 
nosity function. When combined with SDSS data, a 
least-squares fit to a double power-law gives a faint-end 
slope, P = —1.62 ± 0.19, and location of the break at 
M* = -25.6. 

Our binned QLF assumes that all of the QSO can- 
didates are at z = 3.2, which may skew the fitted pa- 
rameters because of luminosity functi on evolution over 
our redshift range and Eddington Bias (jEddingtonll 1 9 1 3| ) 
due to large dispersions in the actual absolute magnitude 
distribution. Therefore, we have performed a maximum 
likelihood fit of the apparent magnitude distribution of 
our sample with that inferred from a specific QLF over 
this redshift range. Our results are slightly different, with 
a shallower faint-end slope, /3 = —1.43 ± 0.15, and a 
somewhat fainter break at M* ~ —24.9. This fit is more 
accurate as it does not assume that all of the QSOs are 
at the same redshift. 

The fitted slope is consistent, within the errors, with 
values measu red at low redshift (0.5 < z < 2.0), /3 = 
— 1.45 (Rich ards et al.l 12005) and therefore does not re- 
quire evolution in the faint end slope of the luminos- 
ity function. Our QLF predicts significantly more faint 
QSOs than sugge sted with initial measurements at z ^ 3 
(|Hunt et al.ll2004D . Although it is difficult to tell wfth our 
limited spectroscopic sample, some o f the difference be - 
tween our faint-end slope and that of Hunt et al.l (|2004l ) 
may be attributed to an increasing number of narrow- 
line, moderately reddened AGN at fainter UV luminosi- 
ties which were excluded from the H04 sample. 

The QLF exhibits some curvature at all magnitudes 
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and, because of this, the parameters for a double power- 
law fit are degenerate. That is, the position of the 
break, (M*,$*), can be fit at different locations along 
the binned QLF, with appropriate changes in the bright 
and faint end slopes (a,/3). This is especially true at 
high redshift, where the difference between the bright 
and faint end slopes appears to decrease. Therefore, one 
must be careful when assigning physical significance to 
the measured values of these parameters when comparing 
to models. 

The QSOs in our sample span the break in the lu- 
minosity function (0.25i* < L < A.OL*) and thus we 
measure the space density of QSOs that comprise the 
majority (55%) of the QSO UV luminosity density at 
these redshifts. When combined with the SDSS sam- 
ple this percentage is more than 70%. Therefore, large 
extrapolations are not required to estimate the effects of 
undetected QSOs. The integrated UV luminosity density 
at z ~ 3.2 is £1450 = 7.3 x lO^" ergs s ~^ Hz"^ h Mpc"'^ 
Using the scaling relation derived by 'Haar dt &: Madaul 
(1996), we infer a maximum HI photoionization rate by 
QSOs, r = 4.5 X lO-i^s-i. This is about 50% of the total 
IGM HI photionization rate at z = 3, requiring compara- 
ble ionizing flux from either starburst galaxies or redder 
AGN that lie outside our color criteria. 
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z=3 QSO Luminosity Function 



TABLEJ 

SWIRE IR Depths from ISurace et al.I (120051 ^ 



Filter 


Central Wavelength (/^m) 


Depth [fiJy, 5cr) 


IRACl 


3.6 


6 


IRAC2 


4.5 


7 


IRAC3 


5.8 


42 


IRAC4 


8.0 


50 


MIPS24 


24.0 


250 



TABLE 2 
WFS Optical Depths 



Filter 


Central Wavelength (A) 


Width (A) 


■niABiVega)'^ 


Depth (Vega, 5a) 


U 


3560 


600 


0.78 


24.3 


9' 


4857 


1400 


-0.09 


25.2 


r' 


6216 


1380 


0.15 


24.5 


i' 


7671 


1535 


0.40 


23.7 


Z 


9100 


1370 


0.54 


22.1 



^ '"AS (Vega) is the Vega to AB conversion factor where mA_B = ^Vega + 
niABiVega). 



TABLE 3 
The broadband optical/IR 
templa te combined w ith the 

ITelfer et al] 

(120021) / TVanden Berk et al.I (120011) 
composite spectra. 



A (/im) 


h 


0.0302 


2.250 


0.0303 


4.086 


0.0304 


6.034 


0.0305 


5.841 


0.0306 


6.000 



Note. — [The complete version of 
this table is in the electronic edition of 
the Journal. The printed edition con- 
tains only a sample.] 



TABLE 4 
Optical Selection Criteria. 



Sample 


Color Criteria (Vega) 




U -g' > 0.33 


^ ~ 3 


g' - r' < 1.0 




1/ - g' > 3.9 X (g' - r') - 2.0 




g' -r' > 1.241 


2 ~ 4 


r' - i' < 1.146 




g' -r' > 2.178 x (r' - i') + 0.09 



TABLE 5 

Optical/IR Photometry of the 100 z ~ 3 QSO candidates. 
Spectroscopic redshifts are given when available. IZ-band 
mac;nitudes in parentheses are 5(7 limits. Two objects with 
spectroscopic confirmation were moved just outside of our color 

CRITERIA (due TO REVISED PHOTOMETRY) AND ARE NO LONGER PART OF OUR 

primary sample. We list these at the end of the table. 



No. 


Name 


RA 


Dec 




U 


9' 


r' 


i' 


Z 


IRACl 


IRAC2 


IRAC3 


IRAC4 


MIPS24 






[Deg] 


[Deg] 


[Vega] 


[Vega] 


[Vega] 


[Vega] 


[Vega] 


[^Jy] 


[^Jy] 


[/"Jy] 


[;^Jy] 


[/"Jy] 


1 


SWIRE_J155907.09+550325.5 


239.77954 


55.05709 




22.08 


21.37 


21.42 


21.17 


20.86 


12 


20 


44 


46 


<250 


2 


SWIRE_J160123. 47+553750.3 


240.34778 


55.63064 




22.75 


22.19 


21.62 


21.14 


21.34 


20 


20 


<42 


<50 


<250 


3 


SWIRE_J160153. 49+542356. 6 


240.47287 


54.39906 




22.15 


21.37 


21.56 


21.01 


20.36 


99 


100 


137 


138 


593 


4 


SWIRE_J160223. 85+552012. 2 


240.59938 


55.33673 




21.92 


21.48 


21.09 


20.59 


19.95 


28 


29 


50 


57 


336 


5 


bWlRE.J 160318. 38+552703. 2 


240.82657 


55.45089 




21.51 


20.95 


20.39 


20.01 


19.63 


56 


83 


133 


225 


827 


6 


bWlRE.J 160345. 78+542337. 2 


240.94077 


54.39367 




20.58 


20.20 


19.68 


19.23 


19.04 


12 


11 


<42 


<50 


<250 


7 


SWIRE_J160422. 50+535454. 9 


241.09373 


53.91525 




22.01 


21.23 


20.60 


20.18 


20.03 


26 


35 


60 


118 


473 


8 


bWlRE.J 160426. 31+553446. 9 


241.10963 


55.57969 




(23.30) 


21.37 


20.91 


20.67 


20.57 


22 


24 


41 


59 


<250 


9 


SWIRE_J160452. 28+542758. 2 


241.21782 


54.46617 


2.98 


22.17 


21.15 


21.30 


21.13 


20.86 


28 


34 


<42 


67 


<250 


10 


bWlKE.J 160520. 52+552704. 6 


241.33548 


55.45128 


3.30 


22.88 


20.51 


20.26 


20.19 


20.13 


48 


50 


66 


73 


461 


11 


SWIRE_,J160528. 14+560635. 6 


241.36725 


56.10989 




22.88 


22.21 


21.78 


21.52 


20.67 


22 


21 


<42 


45 


712 


12 


SWIRE_,J160606. 33+550412. 3 


241.52638 


55.07008 




22.83 


22.25 


21.89 


21.55 


21.13 


16 


19 


<42 


35 


<250 


13 


bWlRE.J 160617.56+541649.5 


241.57318 


54.28042 




23.60 


22.00 


21.21 


20.74 


21.19 


32 


37 


59 


116 


290 


14 


SWIRE_J160621. 18+552532. 7 


241.58827 


55.42574 




23.12 


21.59 


21.11 


20.68 


20.21 


53 


61 


76 


78 


382 


15 


SWIRE_J160637.88+535008.4 


241.65784 


53.83568 


2.943^ 


20.63 


19.79 


19.47 


18.95 


18.75 


109 


169 


312 


739 


2904 


16 


b WlKlij_J160654. 19+554028.4 


241.72578 


55.67455 


2.97 


22.45 


21.72 


21.20 


20.82 


20.63 


21 


29 


<42 


58 


306 


17 


SWIRE_J160656. 10+535633.4 


241.73373 


53.94261 




(23.33) 


21.53 


21.12 


20.66 


20.57 


32 


33 


<42 


68 


276 


18 


SWIRE_J160724.00+533615.2 


241.85002 


53.60422 




(23.25) 


21.77 


21.12 


20.69 


20.48 


15 


17 


<42 


45 


203 


19 


SWIRE. J 160733. 94+554428. 7 


241.89142 


55.74130 




21.24 


20.39 


19.69 


19.04 


18.45 


112 


138 


228 


375 


995 


20 


SWIRE_J160754. 39+533916. 6 


241.97661 


53.65462 


3.01 


23.29 


22.39 


21.93 


21.49 


21.37 


20 


23 


<42 


<50 


202 


21 


SWIRE_J160758.67+543137.8 


241.99446 


54.52716 




21.46 


20.85 


20.59 


20.19 


19.84 


44 


43 


<42 


80 


<250 


22 


bWlKE.J 160824. 08+542003. 7 


242.10033 


54.33435 




22.98 


22.39 


21.98 


21.71 


21.04 


17 


16 


<42 


<50 


<250 


23 


SWIRE_J160850. 55+545800. 5 


242.21063 


54.96680 




(23.34) 


22.60 


21.92 


21.52 


21.45 


11 


11 


<42 


<50 


<250 


24 


SWIRE_J160907.33+543329.8 


242.28056 


54.55827 




22.85 


21.16 


20.57 


20.33 


20.15 


30 


36 


49 


113 


443 


25 


bWlKE.J 160917.26+553638. 2 


242.32191 


55.61062 




22.85 


21.33 


20.96 


20.64 


20.21 


50 


62 


57 


87 


305 


26 


SWIRE_J160947.86+552542.9 


242.44943 


55.42857 




20.47 


19.91 


19.35 


19.03 


18.73 


71 


84 


82 


208 


690 


27 


bWlKE.J 161008. 08+552944. 2 


242.53366 


55.49561 




22.56 


21.73 


21.46 


21.17 


20.92 


14 


14 


<42 


<50 


<250 


28 


SWIRE_J161008. 34+533254. 9 


242.53476 


53.54857 


3.20 


22.81 


21.51 


20.99 


20.18 


19.78 


63 


66 


82 


150 


494 


29 


bWlKE.J 161051. 96+531004. 4 


242.71651 


53.16788 




22.79 


21.82 


21.09 


20.56 


20.15 


32 


37 


58 


110 


417 


30 


bWlRE.J 161 115. 37+534029. 2 


242.81406 


53.67478 




22.89 


22.28 


21.77 


20.75 


20.05 


121 


117 


132 


135 


<250 


31 


bWlRE.J 161 128. 36+553409. 9 


242.86815 


55.56943 




(23.08) 


22.35 


21.70 


21.34 


21.27 


15 


24 


<42 


<50 


<250 


32 


SWIRE_J161 128.47+535809.3 


242.86862 


53.96926 




(23.34) 


22.24 


21.94 


21.46 


21.72 


22 


38 


58 


190 


754 


33 


SWIRE_J161132. 10+542312.0 


242.88374 


54.38667 




22.86 


21.23 


20.58 


20.43 


20.26 


48 


56 


70 


131 


299 


34 


bWlKE_J161142. 40+533104.6 


242.92do7 


53.51794 


3.06 


22.42 


21.35 


21.04 


20.77 


20.51 


37 


43 


57 


86 


302 


35 


bWlK.lij_J161202. 93+532346.9 


243.01221 


53.39636 




(23.12) 


21.23 


20.51 


20.20 


19.97 


26 


29 


<42 


<50 


195 


ou 


o vv iixiij_j J.U J-^oi.yoT^oo^ouo.o 




OO. / oy lu 


3 44 






90 Q'\ 


9n 


90 


90 


22 


^49 


<^ou 




37 


SWIRE_J161300.88+544629.6 


243.25368 


54.77489 




23.27 


21.75 


21.74 


21.09 


21.73 


27 


38 


<42 


76 


166 


38 


SWIRE.J161311.87+542403.7 


243.29944 


54.40102 




22.65 


22.18 


21.83 


21.46 


21.27 


29 


32 


53 


43 


237 


39 


SWIRE.J161326.39+530923.0 


243.35994 


53.15638 


2.83 


20.34 


19.88 


19.47 


19.25 


19.04 


68 


94 


124 


217 


574 


40 


SWIRE_J161341. 27+532956. 4 


243.42195 


53.49900 




22.79 


21.95 


21.24 


21.05 


21.26 


23 


31 


<42 


78 


<250 


41 


SWIRE.J161430.84+555133.3 


243.62848 


55.85926 




22.74 


21.85 


21.53 


21.20 


21.41 


22 


28 


48 


73 


208 


42 


SWIRE_J161433.14+533249.6 


243.63808 


53.54710 




21.04 


20.15 


19.55 


19.36 


19.07 


112 


124 


151 


306 


1068 


43 


SWIRE_,J161442.44+554614.2 


243.67683 


55.77060 




22.86 


21.68 


21.12 


20.76 


20.48 


19 


22 


<42 


56 


<250 


44 


SWIRE_J161446.30+555229.5 


243.69290 


55.87486 




(23.19) 


22.24 


21.53 


20.86 


20.76 


28 


38 


<42 


107 


308 


45 


SWIRE.J161508.88+555514.6 


243.78700 


55.92073 




20.55 


19.71 


19.29 


18.92 


18.78 


55 


75 


122 


301 


924 


46 


SWIRE_J161526.80+555217.4 


243.86165 


55.87150 




22.50 


21.96 


21.33 


20.66 


20.19 


67 


78 


117 


237 


810 


47 


SWIRE_J161530.88+555247.1 


243.87868 


55.87976 




22.91 


22.37 


21.97 


21.49 


20.75 


58 


58 


60 


<50 


<250 


48 


SWIRE.J161549.80+540834.9 


243.95749 


54.14304 




(23.18) 


22.25 


21.84 


21.66 


22.05 


9 


11 


<42 


<50 


<250 


49 


SWIRE_J161626.05+535132.9 


244.10854 


53.85915 




(23.42) 


21.48 


20.82 


20.29 


19.93 


43 


45 


63 


98 


286 



TABLE 5 — Continued 



No. 


Name 


RA 


Dec 








;■' 




Z 


IKACl 


1KAC2 


1RAC3 


IRAC4 


M1PS24 










[\c-ga| 




[\c-ga| 




l\c-ga| 


./'■)y| 










50 


SWIRE_J161634.24+553528.2 


244.14265 


55.59116 




20.14 


19.68 


19.29 


19.01 


18.86 


55 


66 


89 


161 


705 


51 


SWIRE_J161638. 27+555701.4 


244.15947 


55.95039 




21.60 


20.82 


20.42 


19.93 


19.80 


46 


62 


72 


153 


616 


52 


SWIRE_J161704.50+541200.4 


244.26875 


54.20012 




23.24 


22.13 


21.48 


21.15 


21.08 


8 


10 


<42 


<50 


265 


53 


SWIRE_J161719. 00+540154.3 


244.32916 


54.03175 




22.66 


22.05 


21.42 


21.17 


20.98 


28 


40 


45 


107 


356 


54 


sWlRE_J161735. 03+543830.8 


244.39594 


54.64189 




23.66 


22.42 


21.84 


21.64 


21.32 


16 


21 


<42 


44 


381 


55 


SWIRE_J161735. 16+541405.9 


244.39650 


54.23497 




22.68 


20.08 


19.57 


19.22 


19.11 


56 


64 


117 


185 


717 


56 


SWIRE_J161849. 23+543658.3 


244.70511 


54.61620 




22.37 


21.91 


21.53 


21.28 


21.13 


15 


17 


<42 


38 


<250 


57 


SWIRE_J161936. 10+541701. 2 


244.90041 


54.28368 




20.10 


19.74 


20.25 


20.10 


19.84 


195 


235 


329 


378 


1212 


58 


SWIRE_J161959. 51+551453. 8 


244.99796 


55.24828 




(23.29) 


22.05 


21.78 


21.87 


21.48 


12 


12 


<42 


42 


283 


59 


SWIRE_J162004. 22+545023. 4 


245.01758 


54.83984 




22.80 


20.68 


20.45 


20.22 


20.12 


4.5 


56 


104 


224 


989 


60 


SWIRE_J163032. 07+405733. 9 


247.63361 


40.95943 




21.05 


20.70 


20.12 


19.69 


19.49 


11 


10 


<42 


<50 


<250 


61 


SWIRE_J163130. 11+403555. 6 


247.87546 


40.59879 




22.43 


21.43 


21.94 


20.92 


20.29 


113 


99 


117 


110 


650 


62 


SWIRE_J163219. 12+404637.2 


248.07965 


40.77701 




20.26 


19.78 


19.46 


19.14 


19.11 


76 


103 


148 


291 


1032 


63 


SWIRE_J163259. 1 1+401056.6 


248.24628 


40.18239 




22.88 


22.17 


21.76 


21.58 


21.19 


20 


23 


<42 


60 


307 


64 


SWIRE_J163340. 14+404733.0 


248.41725 


40.79249 




22.65 


22.26 


21.99 


21.40 


21.02 


26 


33 


<42 


<50 


<250 


65 


SWIRE_J163343. 54+403739. 5 


248.43143 


40.62764 




20.91 


20.34 


19.71 


19.31 


19.14 


36 


34 


<42 


<50 


<250 


66 


SWIRE_,J163357.79+400225.5 


248.49081 


40.04041 




(23.18) 


21.44 


20.83 


20.56 


20.57 


34 


37 


48 


56 


<250 


67 


SWIRE_J163359. 28+410921.1 


248.49701 


41.15587 




21.52 


21.17 


20.81 


20.42 


19.93 


40 


44 


65 


108 


<250 


68 


SWIRE. J 163403. 26+403845 . 1 


248.51360 


40.64585 




22.81 


22.17 


21.54 


20.93 


20.91 


20 


19 


<42 


<50 


<250 


69 


SWIRE_J163413. 96+412028.3 


248.55818 


41.34119 




22.44 


21.41 


20.92 


20.54 


20.53 


20 


26 


48 


98 


376 


70 


SWIRE_J163417.97+410531.9 


248.57487 


41.09219 




21.12 


20.25 


20.17 


19.60 


19.48 


63 


73 


107 


201 


606 


71 


SWIRE_J163423. 23+400244.0 


248.59680 


40.04555 




20.89 


20.56 


20.04 


19.52 


19.40 


47 


69 


113 


252 


1080 


72 


SWIRE_J163511. 20+404335.4 


248.79665 


40.72651 




21.62 


21.04 


20.64 


19.94 


19.75 


39 


41 


43 


86 


180 


73 


SWIRE_J163527.24+395907.4 


248.86349 


39.98538 




21.89 


21.33 


20.85 


20.45 


20.25 


40 


47 


74 


171 


590 


74 


SWIRE. J 163536.67+412338.7 


248.90279 


41.39407 




22.65 


22.02 


21.77 


21.22 


21.54 


7 


12 


<42 


<50 


<250 


75 


SWiKE.J 163537.37+414904. 2 


248.90570 


A ^ O 1 A 

41.81784 




22.23 


21.00 


20.55 


20.17 


20.14 


29 


36 


63 


126 


415 


76 


SWIRE_J163553. 80+412641. 9 


248.97417 


41.44498 




22.90 


21.96 


21.31 


21.03 


20.60 


25 


32 


53 


87 


<250 


77 


SWIRE_J163604. 98+410307. 7 


249.02077 


41.05215 




22.26 


21.82 


21.21 


20.71 


20.57 


13 


15 


<42 


<50 


<250 


78 


SWIRE_J163627.59+405153.2 


249.11494 


40.86479 




22.54 


21.16 


20.65 


20.18 


19.99 


28 


38 


37 


137 


653 


79 


SWIRE_J163627.66+404218.8 


249.11525 


40.70521 




21.68 


21.20 


20.63 


20.19 


19.98 


14 


14 


<42 


<50 


<250 


80 


SWIRE_J163631. 32+412904. 7 


249.13049 


41.48464 




22.14 


21.71 


21.11 


20.21 


19.78 


208 


313 


332 


451 


1554 


81 


SWIRE_J163657. 12+412850.0 


249.23801 


41.48055 




21.90 


21.55 


21.17 


20.78 


20.33 


17 


26 


<42 


58 


<250 


82 


SWIRE_J16371 1.53+415912. 2 


249.29805 


41.98671 




(23.38) 


22.51 


21.96 


21.59 


21.54 


7 


10 


<42 


<50 


<250 


83 


SWIRE_J163723. 50+414757. 5 


249.34790 


41.79930 




22.86 


21.86 


21.46 


21.09 


21.15 


19 


27 


69 


<50 


<250 


84 


SWIRE_J163733. 22+4131 16. 4 


249.38840 


41.52121 




(23.30) 


21.88 


21.80 


21.15 


21.03 


14 


19 


<42 


37 


<250 


85 


SWIRE_J163744. 75+414245. 6 


249.43646 


41.71266 




21.17 


20.64 


20.10 


19.52 


19.46 


32 


38 


71 


151 


634 


86 


SWIRE_J163822. 19+403650.9 


249.59245 


40.61413 




21.94 


20.97 


20.66 


20.40 


20.40 


28 


28 


<42 


52 


<250 


87 


SWIRE_J163834.42+410014.9 


249.64342 


41.00415 




22.44 


21.75 


21.62 


21.65 


21.29 


13 


14 


35 


52 


264 


88 


SWIRE_J163852. 28+410923.7 


249.71785 


41.15657 




22.43 


22.06 


21.79 


20.89 


20.62 


56 


85 


108 


169 


476 


89 


SWIRE_J163853. 36+421433.2 


249.72234 


42.24255 




20.95 


20.56 


20.17 


20.05 


20.03 


95 


133 


209 


309 


1023 


90 


SW1RE_J163916. 07+414823.7 


249.81694 


41.80657 




(23.27) 


22.06 


21.32 


21.00 


20.61 


18 


23 


<42 


<50 


<250 


91 


SWIRE_J163918. 81+412206.3 


249.82837 


41.36841 




(23.43) 


20.89 


20.36 


19.88 


19.85 


43 


52 


84 


176 


441 


92 


bWiHhv-J 163920.04+421745.0 


249.83350 


42.29605 




22.80 


21.83 


21.31 


21.10 


20.84 


32 


35 


<42 


<50 


<250 


yo 
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4U.OZ (iJiS 




zz.UD 




on '7Ci 




on OK 


0*7 

Z t 


QO 


Oi 


by 


^OKn 

<ZOU 


94 


^WTTJF, Tlfi'^QzLI 7^~[-df\'^QC\Q R 
o vv iJX-Cj—j j.uoy±-L. 1 o^iuoyuy .o 








22.53 


22.19 


21.70 


21.27 


20.98 


g 




<^^z 


<50 


<250 


95 


SWIRE_J163955. 93+410631. 8 


249.98306 


41.10884 




(23.18) 


22.27 


21.80 


21.22 


20.91 


20 


33 


49 


131 


527 


96 


SWIRE.J163956.76+404550.7 


249.98650 


40.76407 




(23.16) 


21.68 


21.31 


99.00 


99.00 


36 


39 


70 


100 


<250 


97 


SWIRE_J164022.80+411548.2 


250.09499 


41.26338 


3.064'* 


20.28 


19.72 


19.20 


18.88 


18.76 


100 


126 


187 


350 


1021 


98 


SWIRE_J164036.64+414916.2 


250.15268 


41.82116 




20.81 


20.46 


20.28 


20.02 


19.69 


39 


51 


99 


128 


380 


99 


SWIRE_J164127.88+412636.4 


250.36617 


41.44345 




22.34 


21.71 


21.26 


21.00 


20.62 


23 


26 


52 


60 


160 


100 


SWIRE.J164229.64+405354.3 


250.62350 


40.89842 




(23.08) 


22.59 


21.98 


21.65 


21.19 


56 


63 


55 


64 


285 


101" 


SWIRE.J160755.85+534020.4 


241.98271 


53.67234 


3.07 


22.46 


22.23 


21.51 


21.27 


21.02 


22 


27 


<42 


54 


185 


102'' 


SWIRE_J161046.25+532540.5 


242.69270 


53.42791 


2.98 


21.41 


21.11 


20.68 


20.51 


20.71 


36 


35 


54 


46 


241 
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TABLE 5 — Continued ^ 
O 



No. 


Name 


RA 


Dec 




u 




;■' 




z 


IRACl 


1KAC2 


1RAC.3 


IRAC4 MIPS24 












[\c-ga| 




[\c-ga| 




[\c-ga| 








[/'•ly. l/'-iy, 



Note. — Typical uncertainties are ~ 0.04 mags in the optical and 10% in infrared fluxes. 

a b 

Zspec from SDSS. Revised photometry moved object outside of color criteria. 



TABLE 6 

Redshifts and Optical/IR Photometry of the seven 
spectroscopically confirmed ^ ~ 4 (g'-dropoux) qsos and two sdss 
qsos in our sample. [/-band magnitudes in parentheses are 5cr limits. 



No. 


Name 


RA 


Dec 


^spec 


U 


g' 


r' 


i' 


Z 


IRACl 


IRAC2 


IRAC3 


IRAC4 


MIPS24 






[Dcg] 


[Dcg] 


[Vega] 


[Vega] 


[Vega] 


[Vega] 


[Vega] 






[^Jy] 




[^Jy] 


1 


oWlRiij_JiDUooz.io+o4zooi.U 


z4i.oooy4 


04.44iy4 


3.82 


23.13 


22.24 


20.90 


20.74 


20.25 


32 


36 


65 


106 


387 


2 


SWIRE.J160705.15-I-533558.7 


241.77147 


53.59965 


3.653'^ 


(22.82) 


19.47 


18.18 


17.89 


17.55 


374 


412 


665 


1392 


5727 


3 


SWIRE_J160907.51-|-535028.0 


242.28130 


53.84110 


3.62 


(22.73) 


22.75 


21.19 


20.69 


20.33 


44 


42 


48 


<50 


231 


4 


S WIRE.J160934. 1 1+550015. 2 


242.39211 


55.00421 


3.52 


22.44 


21.30 


19.79 


20.18 


19.85 


39 


40 


40 


60 


<250 


5 


SWIRE.J161143.22+553157.5 


242.93008 


55.53263 


3.58 


(23.04) 


21.30 


19.99 


19.72 


19.43 


73 


81 


144 


278 


1368 


6 


SWIRE.J161243.17+535827.4 


243.17989 


53.97428 


3.83 


(23.34) 


22.49 


21.02 


20.80 


20.73 


22 


21 


36 


<50 


<250 


7 


SWIRE_J161557.43+545915.8 


243.98929 


54.98772 


3.880 


(23.29) 


23.96 


21.77 


21.20 


21.34 


36 


30 


<42 


55 


<250 


8 


SWIRE.J164326.24+410343.4 


250.85934 


41.06205 


3.873^ 


(23.11) 


21.49 


20.04 


19.81 


19.80 


46 


39 


39 


57 


289 




SWIRE.J104350.94+583029.3 


160.96227 


58.50813 


3.70 




20.95 


19.43 


18.97 




105 


116 


189 


382 


1945 


10'^ 


SWIRE.J160543.05+535829.2 


241.42938 


53.97477 


0.354 


(23.53) 


23.22 


21.48 


20.70 


20.12 


27 * 21 


< 42 


< 50 


<250 






SWIRE.J161159.69+544211.3 


242.99872 


54.70313 


0.390 


22.86 


22.55 


20.72 


20.55 


19.99 


57 


57 


<42 


151 


<250 


I2hd 


SWIRE.J105002.88+574720.0 


162.51199 


57.78890 


? 




21.55 


20.05 


19.73 




215 


391 


677 


1220 


3519 



Note. — Typical uncertainties are ~ 0.04 mags in the optical and ~ 10% in infrared fluxes. ^ 
^ Zspec from SDSS.^ Selected in the SWIRE Lockman field. ^ Low redshift interlopers.^ Could not see emission lines for redshift confirmation. ^ 
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o 
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Siana et al. 



TABLE 7 

Spectral parameter attribute for our QSO SED modeling. 



Parameter 


Mean 


fj 




-0.46 


0.3 


VF(Lya+NV)(A) 


101.7 


25. 



TABLE 8 

Parameters for the four column density ranges used in the Lyq 
Forest simulations. 



Name 


\og{NHi{cm^)) 


No 


K 


7 


b(km/s) 


Lya Forest #1^ 


12-14 


181.36 


1.46 


1.29 


30 


Lya Forest #2=^ 


14-17.2 


1.297 


1.46 


3.10 


30 


Lyman Limit Systems'" 


17.2-20 


0.27 


1.50 


1.55 


70 


Damped Lya Systems'^ 


20-22 


0.055 


1.78 


1.11 


70 



' IKim et ahl ||T997|) '' IStorrie-Lombardi et al] HxM^ " IStorrie-Lombardi fc Wolig l l2000t) 



TABLE 9 
Tabulated Luminosity Function. 



r' M(1450A) 
Vega AB 



n Completeness ^eff 



(10'' Mpc^) (10'^ mag-l Mpc'^) 



18.25 


-27.11 





1.00 


3.22 


7.87 


18.75 


-26.61 





1.00 


3.22 


7.87 


19.25 


-26.11 


7 


1.00 


3.22 


7.87 


19.75 


-25.61 


5 


1.00 


3.22 


7.87 


20.25 


-25.11 


12 


0.99 


3.22 


7.83 


20.75 


-24.61 


20 


0.98 


3.21 


7.73 


21.25 


-24.11 


25 


0.86 


3.18 


6.77 


21.75 


-23.61 


31 


0.74 


3.15 


5.82 



< 0.47 

< 0.47 

-1-1.0 
0.7 



1.78 



1.27; 
3.07; 
5.38 
7.38 
10.66 



-1-0.9 
0.5 
-1-1.2 
0.9 
1.5 
1.2 
1.8 
1.5 
-1-2.3 
1.9 



z=3 QSO Luminosity Function 



TABLE 10 

Parameters for double power-law luminosity function. The binned 
qlf assumes that all qsos are at 2 = 3.2, whereas the m.aximlim 

Likelihood values are from fits to the apparent magnitude 
distribution. values in parentheses denote la errors and values in 
square brackets denote fixed parameters. tlie maximlmum-likeliiiood 

fits that include qlf evolution (both pde and ple) are not 
significantly different from the no evolution parameters. listed at 
the bottom are estimates from previous work. 







a 




P 


A/f* /'A■R^ 


St* yL\j mag ivipc ) 


X 




Binned QLF 


















SWIRE only 


-3.53 (2.9) 


-1.66 (0.88) 


-25.8 (2.3) 


3.5 (11.0) 


2.4 


3 


SWIRE only (fixed a) 




-2.85 




-1.29 (0.85) 


-24.8 (1.1) 


10.4 (12.0) 


2.5 


4 


SWIRE+SDSS 


-.3. 


17 (0. 


58) 


-1,98 (0.17) 


-27.1 (0.7) 


0..^7 (0,65) 


9,9 


12 


SWTRE+SDSS (Iixe([ a) 




-2.85 




-1.02 (0.19) 


-2^.(> (()..3) 


1.53 (2.0) 


11.3 


13 


Maximum Likelihood Fit 


















SWIRE only (fixed a, no evol) 




-2.85 




-1.26 (0.21) 


-25.0 (0.30) 


9.0 (3.0) 






(PDE (l + 2)-3) 




-2.85 




-1.22 (0.22) 


-25.0 (0.29) 


9.3 (2.9) 






(PLE) 




-2.85 




-1.25 (0.21) 


-25.0 (0.30) 


9.2 (2.9) 






SWIRE+SDSS (fixed a, no evol) 




-2.85 




-1.43 (0.15) 


-24.9 (0.15) 


8.5 (2.0) 






(PDE (1 + z)-3) 




-2.85 




-1.41 (0.15) 


-24.9 (0.15) 


8.6 (1.9) 






(PLE) 




-2.85 




-1.42 (0.15) 


-24.9 (0.15) 


8.6 (1.8) 






Previous Studies 


















Hunt ct al. (2004) 


-4.56 (0.51) 


-1.24 (0.07) 


-26.7 


2.4 






Pci (1995) 


-3.52 (0.11) 


-1.64 (0.18) 


-25.8 (0.25) 


6.1 (2.5) 






Bongiorno et al. (2007) 




-3.0 




-1.38 


-25.7 


9.8 







Degrees of freedom in fit- 



